uO  ( 
C vl 

>  i 


AD 


USAAVLABS  TECHNICAL  REPORT  66-30 


SINGLE  ROTOR  HELICOPTER  DYNAMICS  FOLLOWING 
POWER  FAILURE  AT  HIGH  SPEEDS 


By 

D.  Cooper 
K.  C.  Hansen 
T.  I.  Kapiita 

June  1966 


CLEAR  NGHOUSE 
FO°  FEDERAL  *  X  ’'rFTC  AND 

TECHNIOAI 
Hard  copy  1  Mic. 

I 


;ation 


Lo6  t/Asr 


/ 


U.  S.  ARMY  AVIATION  MATERIEL  LABORATORIES 
FORT  EUSTIS,  VIRGINIA 

CONTRACT  DA  44-177-AMv.-l95(T) 

UNITED  AIRCRAFT  CORPORATION 
SIKORSKY  AIRCRAFT  DIVISION 
STRATFORD,  CONNECTICUT  D  D  C 

nVrcr'rpnn  PPfFj 

JUL  14  1366  ; 

LTElU 

B 


Distribution  of  this 
document  is  unlimited 


Disclaime  r  s 


The  findings  in  this  report  are  not  to  be  construed  as  an  official  Depart¬ 
ment  of  the  Army  position  unless  so  designated  by  other  authorized 
documents. 

When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a  definitely  related  Government 
procurement  operation,  the  United  States  Government  thereby  incurs  no 
responsibility  nor  any  obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way  supplied  the  said 
drawings,  specifications,  or  other  data  is  not  to  be  regarded  by  impli¬ 
cation  or  otherwise  as  in  any  manner  licensing  the  holder  or  any  other 
person  or  corporation,  or  conveying  any  rights  or  permission,  to  manu¬ 
facture,  use,  or  sell  any  patented  invention  that  may  in  any  way  be 
related  thereto. 

Trade  names  cited  in  this  report  do  not  constitute  an  official  endorsement 
or  approval  of  the  use  of  such  commercial  hardware  or  software. 


Disposition  Instructions 

Destroy  this  report  when  no  longer  needed.  Do  not  return  it  to 
originator. 


Kcrio*, 

sccrio/ 


I  •'*  fy 

ICfSV 

he 

fa  <"«o«nc:o  ,  ! 

/■Mi we* no*  y  ] 

2  ryf 

f*r  .  . I 

^'iino  «/4MU4,il,n  cod 


0,ST-  /  Mt,  *. 


•r«aMc- 


v 


a* 


■%  -rA ' 

w a'/wiWuii  «’ , .  -1 


DEPARTMENT  OF  THE  ARMY 

U.  S.  ARMY  AVIATION  MATERIEL  LABORATORIES 
FORT  EUSTIS.  VIRGINIA  23604 


This  report  has  been  reviewed  by  the  U.  S.  Army  Aviation 
Materiel  Laboratories  and  is  considered  to  be  technically 
sound . 

As  rotary-wing  aircraft  operations  advance  into  higher  speed 
regimes  such  as  200  knots,  the  dynamic  pressure,  proportional 
to  the  velocity  squared,  becomes  increasingly  £mportant.  This 
report  shows  analytically  that  the  vehicle  dynamics  at  such 
speeds  will  be  determined  in  large  measure  by  fuselage  aero¬ 
dynamics  as  a  consequence  of  high  dynamic  pressure.  The 
aerodynamic  forces  on  the  fuselage  will  determine  how  long 
a  pilot  has  to  recover  from  an  emergency  such  as  powerplant 
failure,  if  he  can  recover  at  all.  Recovery  from  powerplant 
failure  is  shown  to  be  far  easier  where  only  one  engine  in 
multiengine  installations  fails  and/or  where  the  forward 
propulsive  force  required  of  the  rotor  is  minimized,  as  by 
use  of  auxiliary  propulsion. 

The  report  is  published  for  the  exchange  of  information  and 
the  stimulation  of  ideas. 


Task  1D121401A14203 
Contract  DA  44-177-AMC-195(T) 
USAAVLABS  Technical  Report  66-30 
June  1966 


SINGLE  ROTOR  HELICOPTER  DYNAMICS  FOLLOWING 
POWER  FAILURE  AT  HIGH  SPEEDS 

SER-50414 


by 

D.  Cooper,  K.  C.  Hansen,  T.  T.  Kaplita 


Prepared  by 


United  Aircraft  Corporation 
Sikorsky  Aircraft  Division 
Stratford,  Connecticut 


for 


U.  S.  ARMY  AVIATION  MATERIEL  LABORATORIES 
FORT  EUSTIS,  VIRGINIA 


Distribution  of  this 
document  is  unlimited 


BLANK  PAGE 


ABSTRACT 


An  exploratory  study  has  been  conducted  to  investigate  the  effects  of  a 
loss  of  power  on  a  single  rotor  helicopter  at  high  speeds.  A  normal  mode 
flexible  blade  analysis  was  combined  with  generalized  helicopter  equations 
of  motion  on  a  hybrid  computer.  Six  mode  shapes  represented  blade  flexi¬ 
bility,  and  nonlinear  differential  equations  of  motion  described  the  heli¬ 
copter  and  blade  motions  and  loading  in  fifteen  degrees  of  freedom.  This 
analytical  method  was  verified  by  correlation  with  flight  test  data  at 
143  knots.  The  method  was  then  used  to  predict  the  effects  of  rotor  de¬ 
sign  variables  and  corrective  control  inputs  on  the  transient  response  of 
a  helicopter  following  a  full  and  partial  loss  of  main  rotor  power  at 
200  and  235  knots.  This  report  presents  the  results  of  that  study. 

A  full  loss  of  power  can  seriously  affect  the  safety  of  a  200  -  235  knot 
helicopter  unless  there  is  a  rapid  and  appropriate  application  of  main 
rotor  and  tail  rotor  controls.  The  power  loss  is  less  critical  when  a 
partial  loss  of  power  occurs  or  for  a  reduced  rotor  propulsive  force. 
Hence,  the  time  available  for  corrective  action  is  increased. 

Rotor  propulsive  force,  hinge  offset,  and  rotor  blade  inertia  were  found 
to  be  important  design  considerations  which  affect  helicopter  and  blade 
motion  and  loading  following  a  loss  of  power.  Blade  pitch-motion  feedback 
concepts  were  not  effective  in  reducing  the  aircraft  motion.  Large  verti¬ 
cal  and  horizontal  tail  surfaces,  however,  provide  improved  stability  in 
the  event  of  an  abrupt  loss  of  power.  The  aircraft  response  is  due  pri¬ 
marily  to  the  fuselage  aerodynamic  characteristics  associated  with  the 
sideslip  and  yaw  motion  generated  by  the  unbalanced  tail  rotor  yawing 
moment.  The  weathercock  stability  provided  by  a  large  vertical  tail  will 
reduce  the  sideslip  angle  and  resultant  rolling  motion  due  to  fuselage 
and  rotor  dihedral  effect.  The  stable  longitudinal  motion  observed  in 
this  study  demonstrates  the  advantages  of  a  large  horizontal  tail. 

The  effects  of  the  rotor  design  variables  on  aircraft  and  blade  motion 
following  a  full  loss  of  power  are  summarized  in  Table  I.  In  the  table, 
the  rotor  systems  are  compared  with  the  articulated  rotor.  The  design 
variables  within  each  box  are  listed  in  descending  order  of  importance. 

Although  these  qualitative  trends  were  established  using  the  S-61F  (NH-3A) 
helicopter  which  has  very  good  longitudinal  and  directional  stability, 
they  are  expected  to  apply,  in  general,  to  less  stable  helicopters. 

Further  studies  are  recommended  to  investigate  additional  design  variables 
and  control  inputs,  and  to  conduct  research  on  the  pilot  delay  time  prob¬ 
lem  to  provide  experimental  data  on  reactions  of  an  unalerted  pilot  to 
power  loss  in  flight. 
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TABLE  I 

EFFECT  OF  ROTOR  DESIGN  VARIABLES  ON 


AIRCRAFT  AND  BLADE  RESPONSE  AT 
200  KNOTS  WITH  A  STABLE  FUSELAGE 


Response 

Characteristic 

1  Effect  of  Rotor  Design  Variable 

Beneficial 

Negligible 

Detrimental 

Rotor  Speed 
Decay 

Low  Lock  Number 

Hingeless  Rotor 
Floating  Hub 

Zero  Offset 

Delta -three 

Tip  Path  Plane 
Feedback 

Fuselage 

Roll 

Amplitude 

Hinge less  Rotor 

Low  Lock  Number 

Floating  Hub 

Zero  Offset 
Delta-three 

Tip  Path  Plane 
Feedback 

Blade 

Flapping 

Increase 

Floating  Hub 

Low  Lock  Number 
Delta-three 

Tip  Path  Plane 
Feedback 

Zero  Offset 

Flapwise 

Vibratory 

Moment 

Hingeless  Rotor 
Low  Lock  Number 
Delta -three 

Tip  Path  Plane 
Feedback 

Zero  Offset 

Chordwise 

Vibratory 

Moment 

Zero  Offset 

Hingeless  Rotor 
Low  Lock  Number 
Delta -three 

Tip  Path  Plane 
Feedback 

Torsion 

Vibratory 

Moment 

Delta -three 

Tip  Path  Plane 
Feedback 
Hingeless  Rotor 
Zero  Offset 

Low  Lock  Number 

Fuselage 

Pitch  j 

Amplitude 

Low  Lock  Number 
Delta -three 

Tip  Path  Plane 
Feedback 

Zero  Offset 
Hingeless  Rotor 

Floating  Hub  ’ 
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FOREWORD 


This  report,  conducted  by  Sikorsky  Aircraft  with  the  cooperation  of  UAC 
Research  Laboratories,  describes  an  analytic  investigation  of  the  effects 
of  power  loss  on  a  helicopter  at  high  speeds.  The  program  was  sponsored 
by  the  U.  S.  Army  Aviation  Materiel  Laboratories,  Fort  Eustis,  Virginia, 
under  the  technical  direction  of  Mr.  R.  Piper. 

Appreciation  is  due  to  Messrs.  E.  S.  Carter,  Jr.,  D.  S.  Jenney,  K.  C.  Mard, 
and  D.  J.  Carnese,  of  Sikorsky  Aircraft,  for  their  guidance.  Particular 
credit  is  due  to  Mr.  P.  Acidiacano  of  the  Aerophysics  Group  of  UAC  labora¬ 
tories  for  the  development  of  the  equations  of  motion.  Messrs.  G.  Paquette, 
E.  A.  Galligan,  D.  H.  Sorenson,  and  Miss  K.  A.  Derek  of  the  Analog  Com¬ 
puting  Group  of  the  UAC  Laboratories,  are  also  commended  for  their  assis¬ 
tance  in  programing  and  conducting  this  study. 
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INTRODUCTION 


Present  day  helicopters,  which  cruise  in  the  100  to  150  knot  speed  regime, 
can  safely  negotiate  a  complete  loss  of  power.  However,  the  larger  per¬ 
turbations  created  by  a  sudden  power  plant  failure  in  future  high  speed 
helicopters  might  result  in  excessive  rotor  stresses,  both  because  of  the 
larger  applied  torque  change  and  because  of  the  decrease  in  rotor  stabil¬ 
ity  associated  with  the  higher  speed.  Although  it  is  reasonable  to  expect 
the  problem  to  be  most  severe  for  helicopters  that  obtain  all  of  their 
propulsion  and  lift  from  the  rotor  system,  it  might  also  exist  for  high 
speed  helicopters  that  derive  some  propulsive  force  from  auxiliary  power 
plants . 

In  order  to  reduce  the  conjecture  regarding  the  seriousness  of  the  stress 
and  blade  motion  following  an  abrupt  power  plant  failure  on  a  high  speed 
helicopter,  Sikorsky  Aircraft  has  conducted  an  exploratory  analytical 
investigation  under  contract  to  the  U.  S.  Army  Aviation  Materiel  Labora¬ 
tories  (Contract  DA  44-177-AMC-195(T) ) .  The  goal  of  this  study  was  to 
determine  the  influence  of  various  combinations  of  design  parameters, 
flight  conditions,  and  corrective  control  input  patterns  on  the  blade 
stress  and  motion  of  a  helicopter  following  complete  and  partial  power 
loss  at  high  speeds.  This  report  gives  results  of  that  study. 

The  prediction  of  helicopter  and  rotor  dynamic  behavior  after  a  power 
plant  failure  involves  consideration  of  a  large  number  of  coupled,  complex, 
and  nonlinear  effects.  The  classical  approach  to  the  analysis  of  heli¬ 
copter  maneuvering  flight  considers  linear  solutions  only.  This  linear 
approach  to  certain  problems  is  still  very  attractive.  These  solutions, 
which  are  useful  for  generating  trends  for  systems  with  feedback  circuits, 
are  not  particularly  suited  to  calculation  of  helicopter  motion  following 
engine  malfunction. 

A  nonlinear  method  for  calculating  the  applied  and  inertia  loads  acting  on 
a  helicopter  during  various  maneuvers  was  developed  employing  actual  force 
and  moment  equations.  This  method  describes  the  aircraft  in  six  degrees 
of  freedom  in  the  body  axis  system  and  the  rotor  in  flapping  and  rotor 
speed  degrees  of  freedom.  With  the  inclusion  of  rotor  blade  inplane  motion 
(lead-lag)  and  intricate  rotor  blade  aeroelastic  characteristics,  the 
method  is  well  suited  to  the  analysis  of  dynamic  response  of  a  helicopter 
following  power  plant  failure  at  high  speeds. 

The  Sikorsky  Normal  Mode  Flexible  Blade  Analysis  system  was  used  for  the 
rotor  system.  It  solves  the  blade  equations  by  time-wise  step-by-step 
integration  and  is  directly  applicable  to  the  calculation  of  rotor  tran¬ 
sient  motion.  This  analysis  represents  the  blade  deflections  in  a  conven¬ 
tional  manner  by  a  summation  of  rvormal  vibratory  blade  elastic  modes  to 
give  the  fully  coupled  flapwise,  chordwise,  and  torsional  bending  response. 
This  analysis  provides  a:a  accurate  representation  of  the  rotor  blade  aero¬ 
elastic  characteristics  and  includes  blade  inplane  (lead-lag)  degree  of 
freedom.  Essential  agreement  with  test  data  has  been  obtained  in 
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correlation  studies. 


The  above  analysis  was  coupled  on  the  hybrid  digital  analog  computer  with 
the  general  helicopter  simulation  force  and  moment  equations  of  motion. 

This  coupling  provides  a  complete  and  accurate  representation  of  body  and 
blade  motions  as  well  as  stress  levels  at  critical  blade  spanwise  stations. 
This  combined  analysis,  representing  a  significant  advancement  in  analyti¬ 
cal  techniques,  is  general  and  well  suited  to  a  variety  of  VTOL  aircraft 
and  maneuvering  flight  conditions.  Good  correlation  has  been  achieved 
with  S-6l  flight  test  data. 

This  report  presents  the  results  of  the  application  of  this  method  to  the 
analysis  of  helicopter  rotor  and  fuselage  motion  following  an  exponential 
reduction  of  50$  and  100$  in  power  supplied  to  the  rotor  system  and  subse¬ 
quent  insertions  of  corrective  rotor  control  inputs.  The  major  portion  of 
this  study  is  concerned  with  flight  speeds  at  200  and  235  knots  of  the 
six-bladed,  S-61F  (NH-3A)  high  speed  research  helicopter,  configured  with 
various  combinations  of  flapping  hinge  offset,  twist,  tip  path  plane  and 
flapping  feedback,  Lock  number,  and  drag-to-lift  ratios. 
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SCOPE  OF  PROGRAM 


The  effects  of  power  loss  on  helicopter  behavior  and  blade  structural  load¬ 
ing  have  been  studied,  for  general  rotor  design  parameters,  at  200  and  235 
knots.  The  influence  of  flapping  hinge  offset,  articulation  method,  twist, 
Lock  number  and  pitch-blade  motion  coupling  has  be  i  investigated  for  50$ 
and  100$  reduction  in  the  available  rotor  system  dr  ving  torque.  In  this 
study,  a  partial  power  loss  or  50$  loss  in  available  rotor  system  driving 
torque  is  considered  as  a  loss  of  one  engine.  A  full  power  loss  assumes 
the  loss  of  both  engines.  Auxiliary  propulsive  force  is  assumed  to  be 
provided  by  an  independent  power  source.  Four  values  of  drag-to-li  ft 
ratios,  D/L,  have  been  considered  in  order  to  determine  the  influence  of 
auxiliary  propulsion.  In  addition,  various  types  of  corrective  control  in¬ 
puts  were  investigated,  with  delay  times  ranging  from  0.2  second  to  0.6 
second. 

In  order  to  establish  qualitatively  the  interrelation  among  design  combina¬ 
tions,  flight  conditions,  and  input  patterns,  hybrid  computer  time  history 
data  of  body  and  blade  motions  were  recorded  for  each  basic  variation.  The 
flight  conditions  investigated  included  the  following  combinations: 

200  knots,  D/L  =  0.15,  =  -4° 

200  knots,  D/L  =  0.075;  #1  =  -i+° 

200  knots,  D/L  =  0,  81  =  0° 

235  knots,  D/L  =  -O.O63,  6±  =  0° 

where  D/L  is  the  ratio  of  rotor  propulsive  force  to  rotor  lift  (drag-to- 
lift  ratio) . 

The  design  combinations  were  incorporated  into  seven  rotor  systems  in¬ 
stalled  on  the  same  base  aircraft.  The  high  speed  S-6lF  research  aircraft, 
described  in  Appendix  I,  was  used  as  a  basis  for  this  study.  A  gross 
weight  of  11,900  pounds  was  assumed. 

The  control  input  portion  of  this  study  was  conducted  for  six  rotor  con¬ 
figurations  and  three  D/L  conditions  to  determine  the  influence  of  rate 
and  magnitude  of  inputs,  and  coupling  of  controls. 

Blade  structural  loading  and  motion  time  history  data  were  recorded  with 
and  without  control  inputs  for  those  cases  considered  to  be  most  severe. 
Chordwise,  flapwise,  and  torsion  bending  moments  have  been  computed  at 
three  radial  stations,  utilizing  realistic  blade  structures.  In  addition, 
chordwise  bending  moments  have  also  been  calculated  in  the  vicinity  of  the 
lag  damper.  The  data  have  been  analyzed  to  define  qualitatively  the  ef¬ 
fects  of  the  design  parameters  and  control  inputs  on  blade  loading. 
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DESCRIPTION  OF  ROTOR  SYSTEMS 


The  basic  design  combinations  of  flapping  hinge  offset,  articulation  method, 
twist.  Lock  number,  and  pitch-blade  motion  coupling  were  investigated  by 
means  of  the  following  seven  rotor  systems: 

a.  Fully  articulated  six-bladed  rotor  with  a  flapping  hinge  offset 
ratio  (e/R)  of  0.0339,  X*  10.62. 

b.  Rotor  system  (a)  above  with  individual  blade  pitch-flap  coupling 
delta-three  =  30°  • 

c.  Rotor  system  (a)  above  with  individual  blade  pitch  coupling  with 
the  integrated  flapping  of  all  blades. 

d.  Rotor  system  (a)  above  with  a  flapping  hinge  offset  ratio 
(e/R)  of  0. 

e.  Rotor  system  (a)  above  with  a  decreased  Lock  number, 

X=  7-88. 

f .  Rotor  with  four  low  aspect  ratio  blades  rigidly  attached  to 
the  hub. 

g.  Rotor  system  (f)  above  with  a  floating  hub. 

The  value  of  delta-three,  used  for  this  study  for  rotor  b,  is  a  typical 
value  considered  for  use  on  present  day  high  speed  helicopters  to  reduce 
the  gust  sensitivity  of  the  main  rotor .  A  comparable  amount  of  tip  path 
tilt  feedback  was  used  for  rotor  c.  For  rotor  c,  the  tilt  of  the  tip  path 
plane  was  sensed  by  feeding  back  blade  flapping  velocity  to  blade  pitch.  A 
stable  feedback  was  used  for  both  rotors. 

The  physical  properties  of  the  various  rotor  systems  are  given  in  Appendix 

I. 

The  study  conditions  investigated  for  each  of  those  rotor  systems  are  out¬ 
lined  in  Table  II,  and  the  rotor  operating  conditions  are  listed  in  Table 
III.  With  the  advancing  blade  t'p  Mach  number  specified  contractually,  as 
shown  in  Table  III,  a  maximum  rotor  tip  speed  of  633  fps  was  selected  as 
providing  the  best  performance  characteristics.  The  most  critical  opera¬ 
ting  condition,  in  terms  of  stall,  is  at  a  drag-to-lift  ratio  (d/l)  of  0.15 
at  200  knots.  This  condition  sized  each  rotor  so  that  a  maximum  value  of 
rotor  lift  coefficient  (Cl/ct)  of  0.04A  was  available  within  the  "upper 
stall"  limit  of  (bCqd/^)  =  0.008.  This  is  a  reasonable  limit  for  steady 
flight  operation.  For  the  fully  articulated  rotor  with  six  blades,  this 
operating  condition  dictated  a  rotor  lift  of  11,900  pounds.  At  235  knots, 
however,  a  D/L  =  -  O.063  was  selected,  in  order  to  provide  a  reasonable 
amount  of  rotor  torque  within  an  "onset  of  stall"  limit  of  (bCq^ ;(r)  =  0.00^+ . 
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TABLE  II 
STUDY  CONDITIONS 


Forward.  Speed 

200  Knots 

235  Knots 

Blade  Twist 

-4  Degrees 

0  Degrees 

Rotor  System 

0.15 

0.075 

0. 

-0 . 063 

a  Fully  Articulated 

X 

X 

X 

X 

b  Delta -3  Feedback 

X 

- 

X 

X 

c  Tip  Path  Plane  Tilt 
Feedback 

X 

_ 

X 

d  Zero  Offset 

X 

- 

X 

- 

e  Low  Lock  Number 

X 

X 

X 

X 

f  Hingeless  Rotor 

X 

X 

X 

X 

g  Floating  Hub 

X 

X 

X 

- 

TABLE  III 

ROTOR  OPERATING  CONDITIONS 


V 

mt90° 

floR 

H- 

(CL/  *  )stall  (D/L) 

Knots 

- 

fps 

- 

- 

- 

200 

0 .870 

633 

0.534 

0.044 

0.15 

200 

0 .870 

633 

0.534 

0.052 

0.075 

200 

0.870 

633 

0.534 

0.058 

0. 

235 

0.923 

633 

0.627 

0.044 

-O.O63 

vi1 

'V 

— i 

i 

« 
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DEVELOPMENT  OF  COMPUTER  PROGRAM 


The  prediction  of  helicopter  and  rotor  system  dynamic  response  following  a 
power  plant  failure  involves  consideration  of  a  large  number  of  coupled, 
complex,  and  nonlinear  effects.  The  analysis  used,  programed  on  a  hybrid 
computer  system,  describes  the  aircraft  in  six  degrees  of  freedom  in  the 
body  axis  system  and  the  rotor  in  flapping,  lagging,  bending,  and  rotor 
speed  degrees  of  freedom,  and  represents  the  blades  by  a  summation  of  nor¬ 
mal  vibratory  blade  elastic  modes  to  give  a  fully  coupled  flapwise,  chord- 
wise,  and  torsional  bending  response. 

The  program  includes  the  following: 

1.  Six  degree  of  freedom  fuselage  equations  of  motion  with  no  small  angle 

assumptions . 

2.  Nonlinear  fuselage  aerodynamic  wind  tunnel  data. 

3*  Nonlinear  rotor  force  and  moment  equations  of  motion  for  each  blade 
with  no  angle  of  attack  or  advance  ratio  restrictions. 

4.  Nonlinear  rotor  blade  airfoil  section  aerodynamic  data  including  stall 
and  Mach  number  effects. 

5-  Individual  blade  flapping  and  lagging  degrees  of  freedom. 

6.  Rotor  speed  degree  of  freedom. 

7.  Six  aerodynamically  coupled  blade  deflection  degrees  of  freedom. 

A.  Three  Flapwise  Modes 

B.  Two  Chordwise  Modes 

C.  One  Torsional  Mode 

In  the  analysis,  tail  rotor  power  is  assumed  small  and  is  not  included  in 
the  decay  of  main  rotor  speed.  For  the  aircraft  trimmed  for  zero  roll  in 
the  200- knot  range  of  forward  speed,  tail  rotor  power  is  only  about  3 % 
of  the  main  rotor  power  and,  therefore,  is  a  second  order  effect. 

In  the  blade  aeroelastic  part  of  the  program,  although  one  reference  blade 
was  monitored,  a  complete  aerodynamic  and  aeroelastic  solution  was  obtained 
for  each  blade.  This  procedure  insured  accurate  calculation  of  rotor  per¬ 
formance  and  simplified  the  analysis  of  the  floating  hub  rotor  system. 

A  description  of  equations  of  motion,  flexible  blade  simulation,  assump¬ 
tions  and  method  of  solution  is  given  in  Appendix  II.  A  list  of  equations 
used  in  this  study  is  given  in  Appendix  III. 
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PROGRAM  CHECK  OUT  AND  CORREIATION 


Verification  of  the  hybrid  computer  program  was  conducted  by  correlating 
calculated  data  with  actual  flight  test  data  of  a  simulated  power  failure. 
Flight  test  data  of  the  Sikorsky  CH-3C  aircraft,  described  in  Appendix  I, 
was  utilized  in  this  study.  A  comparison  with  flight  test  data  of  the  cal¬ 
culated  response  of  the  aircraft  due  to  a  simulated  power  loss  at  143  knots 
is  presented  in  Figure  1.  Actual  pilot  control  inputs,  shown  in  Figure  2, 
as  well  as  an  exponential  decay  of  power  with  time  w^re  used  in  the  com¬ 
puter  program.  Engine  power  was  permitted  to  decay  exponentially  to  35$  of 
its  trim  value  in  0.5  second.  However,  as  seen  in  Figure  1,  after  1.4  sec¬ 
onds,  the  power  lever  was  changed  to  provide  more  power.  This  power  rise 
was  not  simulated  in  the  computer  program  due  to  the  limited  information  on 
the  engine  response  characteristics  and  the  actual  power  change. 

In  Figure  1,  the  calculated  response  of  the  aircraft  agrees  favorably  with 
the  flight  test  data  during  the  time  interval  that  the  power  lever  was  held 
fixed.  In  particular,  calculated  rotor  speed  and  blade  lag  angle  show 
excellent  correlation  with  flight  test  data.  Fuselage  angular  rates  and 
displacements  show  good  correlation  in  both  trends  and  magnitudes. 

In  previous  correlation  studies  conducted  in  development  of  stability  and 
control  programs,  a  significant  variation  in  the  response  of  unstabilized 
helicopters  was  observed  whenever  the  initial  conditions  were  varied  or 
whenever  mild  atmospheric  turbulence  was  introduced.  It  is  reasonable  to 
assume  that  during  this  test  the  aircraft  was  not  in  as  perfect  trim,  nor 
in  as  ideal  weather  conditions  as  can  be  commanded  on  the  computer.  These 
non-ideal  conditions  will  affect  the  aircraft  response  following  the  power 
reduction.  The  long  term  difference  between  calculated  and  test  value  of 
the  aircraft  angular  rates  and  displacements  is,  therefore,  attributed  to 
the  difference  in  these  conditions. 
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Figure  1.  Comparison  of  Calculated  Aircraft  Response  Due  to  a  Simulated 
Power  Loss  With  Flight  Test  Data  at  1^3  Knots 
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DISCUSSION  OF  RESULTS 


Time  history  data  of  body  and  blade  motions  were  recorded  for  each  rotor 
configuration  in  order  to  establish  the  interrelation  among  rotor  paramet¬ 
ers,  flight  conditions,  and  corrective  control  input  patterns.  A  total  of 
178  combinations  were  investigated  and  the  most  significant  have  been  ana¬ 
lyzed  and  discussed  here.  A  list  of  the  configurations  and  flight  condi¬ 
tions  is  given  in  Table  II,  and  a  selected  list  of  the  corrective  control 
inputs  discussed  in  this  section  is  presented  in  Table  IV.  As  seen  in  the 
latter  table,  these  inputs  include  various  combinations  of  step  and  ramp 
pedal  inputs  and  ramp  inputs  of  cyclic  and  collective  blade  pitch.  All 
ramp  inputs  are  applied  at  the  rate  of  ^ .  5  degrees  per  second.  The  trim 
condition  for  each  rotor  system  is  given  in  Table  V. 

Typical  time  history  data  recorded  for  this  study  are  presented  in  Figures 
3  and  4.  The  response  of  a  fully  articulated  rotor,  rotor  (a),  following 
a  full  power  loss  at  200  knots,  (D/L  =  0. 15)  with  no  control  input  is  pre¬ 
sented  in  Figure  3.  The  effect  of  a  full  corrective  control  input,  6b, 
initiated  0.2  second  after  the  power  loss  is  shown  in  Figure  4. 

In  this  section,  helicopter  and  blade  motion  following  a  loss  of  power  is 
discussed  and  then  the  blade  loading  associated  with  these  disturbances  is 
considered.  The  effects  of  Lock  number,  flapping  hinge  offset,  and  blade 
pitch-blade  motion  feedback  on  rotor  speed,  flapping,  lagging,  and  fuselage 
roll  angle  with  no  control  inputs  are  treated  first.  The  discussion  is 
based  on  values  recorded  2  seconds  after  power  loss  or  on  values  of  the 
first  peak  after  power  loss,  section  A.  This  analysis,  based  on  discrete 
values,  is  further  supplemented  by  a  detailed  analysis  of  helicopter  and 
blade  transient  response,  section  B.  The  time  history  data  were  replotted 
to  show  the  effects  of  the  power  loss  by  a  relative  comparison  of  rotor 
configurations,  flight  conditions,  and  corrective  control  inputs.  The 
blade  loading  discussion,  section  C,  includes  both  an  analysis  of  the  trim 
bending  moments  for  all  configurations  except  the  floating  hub  rotor  and  an 
analysis  of  the  maximum  vibratory  moment  characteristics  for  several  rotor 
configurations,  control  inputs,  and  main  rotor  drag  to  lift  ratios,  D/L. 

Note  that  D/L  represents  the  ratio  of  rotor  propulsive  force  to  rotor  lift. 
This  section  is  then  canpleted  by  a  discussion  of  the  time  history  of  peak 
and  cne-half  peak-to-peak  bending  moment  envelope  for  various  configurations 
and  control  inputs. 

A.  EFFECT  OF  ROTOR  DESIGN  PARAMETERS  ON  HELICOPTER  AND 
BLADE  MOTION  PEAK  VALUE  CHARACTERISTICS  FOLLOWING 
POWER  LOSS  -  NO  CORRECTIVE  CONTROL  APPLIED 


Effect  of  Lock  Number 


A  2bf0  reduction  in  Lock  number  from  10.62  to  7*88  does  not  significantly 
affect  the  rotor  speed  following  loss  in  power  as  shown  in  Figure  5.  The 
rotor  speed,  as  a  percentage  of  trim  rotor  speed  measured  two  seconds  after 
the  loss  in  power,  is  used  as  a  criterion  for  comparison.  The  rotor  speed 
for  each  configuration  differs  only  by  about  2 f>.  Following  a  full  powei 
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loss,  the  rotor  with  the  low  Lock  number  blades  maintains  a  higher  rotor 
speed  than  the  standard  rotor.  Both  rotors  maintain  cbout  the  same  speed 
following  a  partial  loss  in  power. 

The  higher  inertia  of  the  low  Lock  number  blades  reduces  the  flapping  dis¬ 
turbance  as  seen  in  Figure  6.  The  change,  from  trim,  in  the  one-half  peak- 
to-peak  maximum  vibratory  flapping  amplitude  is  used  for  comparison. 

Neither  rotor  experiences  a  change  in  blade  flapping  following  a  partial 
loss  in  power  for  drag  to  lift  ratios  of  0.075>  0,  and  -O.O63.  Although 
the  high  inertia  blades  would  cause  the  rotor  disk  to  lag  behind  the  air¬ 
craft  roll  motion  and  would  indicate  a  larger  flapping  motion  relative  to 
the  shaft,  the  reduction  in  roll  motion  reduced  the  rotor  disk  lag,  there¬ 
by  reducing  blade  flapping. 

The  change  in  steady  lag  angle  is  also  reduced  as  the  blade  Lock  number  is 
decreased,  as  shown  in  Figure  7-  Tb^  tip  mass  in  the  centrifugal  force 
field  acts  as  a  restoring  moment  to  educe  the  initial  lagging  position  of 
the  blades  and  the  blade  lead- lag  excursion. 

The  increased  hub  moment  damping  provided  by  the  higher  inertia  of  the  low 
Lock  number  blades  reduces  the  aircraft  roll  response,  as  seen  in  Figure  8. 
A  comparison  of  the  first  peak  in  fuselage  roll  angle  after  the  power  loss 
is  presented  in  this  figure.  For  both  configurations,  the  severity  of  the 
roll  disturbance  is  considerably  reduced  as  the  rotor  propulsive  force  re¬ 
quirement  is  reduced,  with  lower  D/L.  At  235  knots,  D/L  =  -O.063,  the  air¬ 
craft  displays  little  or  no  roll  disturbance  even  for  a  full  loss  in  power. 
Therefore,  increasing  the  rotor  blade  moment  of  inertia  will  effectively 
reduce  the  blade  and  aircraft  disturbance  following  a  loss  of  power. 

Effect  of  Blade  Pitch-Blade  Motion  Feedback 


Delta -three  and  tip  path  plane  feedback  will  not  be  beneficial,  as  seen  in 
Figures  9  through  12.  Compared  to  a  standard  fully  articulated  rotor, 
rotor  speed  falls  off  faster,  as  seen  in  Figure  9-  Within  the  first  two 
seconds  following  loss  of  power,  the  delta-three  rotor  loses  more  speed 
than  the  other  configurations.  Flapping  disturbance,  however,  is  decreased 
when  a  full  power  loss  occurs,  and  is  slightly  increased  following  a  part 
power  loss  (Figure  10).  Blade  lead-lag  motion  is  also  higher  with  blade 
motion  feedback  as  seen  in  Figure  11.  The  time  history  data  for  the  delta- 
three  rotor  were  terminated  about  1.5  seconds  after  the  full  power  loss 
(Figure  2l).  The  termination  occurred  because  the  aircraft  motion  exceeded 
reasonable  value  limits  programed  in  the  computer.  Blade  flapping  motion 
was  therefore  not  fully  developed  for  this  flight  condition.  Of  significant 
interest,  however,  are  the  peak  roll  response  characteristics  shown  in 
Figure  12. 

Delta-three  and  tip  pach  plane  tilt  feedback  aggravate  the  roll  transient 
response  of  a  helicopter  following  a  loss  of  power.  Although  the  peak  roll 
amplitudes  are  nearly  the  same  for  a  part  pc  wer  loss,  the  delta-three  rotor 
experiences  a  larger  roll  disturbance  than  either  the  tip  path  plane  feed¬ 
back  or  standard  rotor  following  a  full  loss  of  power.  Both  feedback  rotors 
experience  a  larger  roll  disturbance  than  th^  standard  rotor.  As  rotor 
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propulsive  force  requirements  are  reduced,  with  lower  D/L,  all  rotors  have 
the  same  reduced  response. 

Effect  of  Flapping  Hinge  Offset 

A  reduction  in  flapping  hinge  offset  ratio,  e/R,  with  the  associated  reduc¬ 
tion  in  hub  moment  damping,  also  aggravates  the  transient  response  of  a 
helicopter  following  a  power  loss  at  high  speeds.  This  is  evident  in  Fig¬ 
ures  13  through  16.  Compared  with  a  rotor  with  a  3-4$  flapping  hinge  off¬ 
set  ratio,  a  rotor  with  zero  offset  will  experience  a  faster  decay  in  rotor 
speed  (Figure  13),  a  large  flapping  disturbance  (Figure  14),  and  a  larger 
excursion  in  lead-lag  angle  following  both  full  and  part  power  losses  (Fig¬ 
ure  15)*  The  reduction  in  hub  damping  with  zero  offset  also  causes  the  roll 
disturbances  to  be  larger,  as  seen  in  Figure  16.  The  disturbance,  however, 
is  reduced  for  both  configurations  as  the  rotor  propulsive  force  require¬ 
ments  are  reduced. 

An  analysis  of  the  recorded  time  history  data  indicates  that  Lock  number, 
blade  motion  feedback,  and  flapping  hinge  offset  do  not  significantly  affect 
the  pitch,  yaw,  and  sideslip  motion.  The  peak  values  of  these  quantities 
are  therefore  presented  in  Figure  17  as  a  single  curve.  It  is  seen  in  this 
figure  that  the  pitch,  yaw,  and  sideslip  response  are  primarily  affected  by 
rotor  propulsive  force  requirements,  D/L,  and  fuselage  aerodynamic  charac¬ 
teristics,  as  will  be  shown  in  the  next  section. 

B.  EFFECT  OF  ROTOR  DESIGN  PARAMETERS,  ROTOR  PROPUISIVE 
FORCE,  AND  CORRECTIVE  CONTROL  INPUT'S  ON  HELICOPTER 
AND  BLADE  MOTION  TIME  HISTORY  FOLLOWING  POWER  LOSS 


Although  the  evaluation  of  discrete  values  of  helicopter  and  blade  motion 
presented  above  does  show  some  distinctive  trends,  a  complete  time  history 
of  the  transient  response  is  required  for  a  thorough  appreciation  of  the 
effects  of  rotor  parameters  and  control  inputs.  For  this  reason,  selected 
time  history  data  of  helicopter  and  blade  motions  have  been  replotted  and 
presented  in  Figures  18  through  47.  The  transient  response  data  are 
presented  for  various  combinations  of  rotor  configurations,  corrective  con¬ 
trol  inputs,  and  drag-to-lift  ratios.  The  data  presented  include  blade 
flapping  and  lagging  motion,  rotor  speed,  fuselage  roll  and  sideslip  angles, 
and  fuselage  yaw  and  pitch  rates. 

Part  Power  Loss  -  No  Control  Input 

The  transient  helicopter  and  blade  motion  following  a  part  power  loss  at 
a  D/L  =  0.15,  200  knots,  for  all  rotor  configurations  is  presented  in  Fig¬ 
ures  18  and  19.  The  response  of  a  helicopter  with  a  blade  pitch-blade  flap 
feedback  rotor  system  (rotor  b);  a  blade  pitch-tip  path  plane  feedback 
rotor  system  (rotor  c),  and  a  reduced  Lock  number  rotor  system  (rotor  e)  is 
compared  with  the  response  of  a  fully  articulated  rotor  system  (rotor  a)  in 
Figure  18.  The  response  of  rotor  systems  with  zero  flapping  hinge  offset 
(rotor  d),  hingeless  blades  (rotor  f),  and  floating  hub  (rotor  g)  is  com¬ 
pared  with  the  response  of  rotor  (a)  in  Figure  19 . 
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Throughout  the  recorded  time  period  in  Figure  l8,  Lock  number  does  not 
appreciably  affect  rotor  speed.  Although  rotor  speed  tends  to  converge  to 
approximately  the  same  value  at  the  end  of  the  period,  the  delta-three  rotor 
has  a  higher  decay  rate  and  the  speed  remains  below  all  others.  About  4 
seconds  after  the  partial  loss  of  power,  rotor  speed  becomes  constant  at 
about  90 °Jo  of  its  trim  value.  At  this  new  rotor  speed,  the  power  required 
by  the  main  rotor  matches  the  power  available  from  the  remaining  single 
engine.  Rotor  speed  can  be  increased  by  reducing  forward  speed  or  entering 
a  descending  flight  condition  by  lowering  collective  pitch. 

The  roll  response,  however,  is  quite  severe.  Rotors  (a),  (b),  and  (c)  roll 
left  more  than  40°  in  2.5  seconds  after  the  loss  of  power.  However,  the 
increased  hub  damping  moment  associated  with  the  low  Lock  number  blades 
limits  its  roll  disturbance  to  about  30°*  The  stabilization  concepts  of 
flapping  and  tip  path  plane  feedback  to  blade  pitch  are  seen  to  be  ineffec¬ 
tive,  because  these  systems  experience  a  larger  roll  motion  than  the  stan¬ 
dard  rotor. 

The  rapid  reduction  in  main  rotor  torque  without  an  accompanying  reduction 
in  tail  rotor  thrust  creates  an  unbalanced  yawing  moment  condition  on  the 
fuselage.  Consequently,  the  aircraft  yaws  to  the  left  about  -5°  to  gener¬ 
ate  a  sideslip  angle  of  approximately  +5°  in  2  seconds.  As  seen  in  Figure 
3,  the  aircraft  is  trimmed  at  a  negative  angle  of  attack  of  -4°.  For  this 
combination  of  fuselage  angle  of  attack  and  sideslip,  the  fuselage  rolling 
moment  is  negative  and  yawing  moment  is  positive  (Figure  78).  The  initial 
yaw  motion  to  the  left  is  reversed  in  about  1-5  seconds  by  the  fuselage 
yawing  moment,  and  the  aircraft  rolls  to  the  left  due  to  the  stable  fuse¬ 
lage  and  rotor  dihedral  effect.  The  pitch  and  yaw  rates,  sideslip  angle, 
and  their  period  of  motion,  however,  are  virtually  the  same  for  all  rotor 
configurations.  The  aircraft  response  is,  therefore,  due  primarily  to  the 
fuselage  aerodynamic  forces  and  moments  associated  with  the  sideslip  and 
yaw  motion  generated  by  the  unbalanced  tail  rotor  yawing  moment. 

The  blade  lead-lag  motion  for  all  configurations  has  a  vibratory  amplitude 
of  about  ±1°.  Blade  flapping  is  relatively  constant  and  then  increases 
markedly  at  about  4  seconds.  The  large  flapping  amplitude  of  the  delta- 
three  rotor  in  comparison  with  the  fully  articulated  rotor  is  evident. 

Higher  harmonics  of  flapping  also  occur  in  the  delta-three  trace.  Since 
fuselage  roll  rate  peaks  at  about  4  seconds  (Figure  3);  the  large  flapping 
amplitude  is  attributed  to  the  angle  of  attack  induced  by  the  roll  rate. 

The  transient  response  of  rotor  speed  for  the  hingeless  rotor  is  initially 
the  same  as  for  the  fully  articulated  rotor  (Figure  19)  •  About  2  seconds 
later,  however,  the  hingeless  rotor  speeds  up.  The  speed-up  is  attributed 
to  the  change  in  rotor  angle  of  attack  distribution  induced  by  the  aircraft 
roll  and  pitch  rates.  This  rotor  maintains  a  zero  roll  attitude  for  about 
1  second,  then  rolls  rapidly  to  -80°  in  about  3  seconds,  and  continues  to 
diverge  with  time.  Although  the  gyroscopic  characteristics  of  the  hinge¬ 
less  rotor  initially  limit  the  effects  of  the  power  loss,  the  unbalanced 
tail  rotor  thrust  generates  a  disturbance  on  the  aircraft,  which  produces 
a  left  roll.  The  gyroscopic  effects  generate  a  nose  down  moment  which 
increases  the  negative  angle  of  attack  of  the  fuselage.  A  larger  negative 
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rolling  moment  is  developed  which  increases  the  left  roll  motion,  causing 
the  aircraft  to  diverge. 

The  effect  of  low  hub  moment  damping  is  evident  in  the  motion  of  the  zero 
offset  rotor.  Aircraft  motion  is  larger  and  rotor  speed  decays  faster. 
The  effect  of  the  aircraft  motion  on  flapping,  as  well  as  the  increase  in 
vibratory  lead-lag  motion  as  rotor  speed  falls  off,  is  evident  in  Figure 

19. 


No  automatic  stabilization  has  been  assumed  in  any  of  the  cases  shown.  The 
divergence  could  probably  be  controlled  by  an  auxiliary  stabilization 
device,  but  the  optimization  of  proper  feedback  was  beyond  the  scope  of 
this  study. 

Full  Power  Loss  -  No  Control  Input 

Effects  of  rotor  propulsive  force,  fully  articulated  rotor 

The  effect  of  a  full  loss  of  power  on  the  transient  motion  of  the  fully 
articulated  rotor  for  various  drag-to-lift  ratios  is  shown  in  Figure  20. 

As  D/L  increases,  the  rate  of  decay  of  rotor  speed  also  increases  for  the 
200  knot  condition  (D/L  =  0,  0.075  and  0.15),  because  of  the  increase  in 
power  loading  with  an  increase  in  D/L.  Since  the  rotor  at  D/L  =  -O.O63 
(235  knots)  has  about  the  same  power  loading  as  at  D/L  =  0  (200  knots),  the 
decay  in  rotor  speed  is  about  the  same. 

As  the  D/L  ratio  increases,  the  fuselage  response  becomes  more  severe. 

This  is  due  to  the  fuselage  aerodynamics.  Larger  slip  angles  occur  and, 
therefore,  larger  roll  angles  due  to  fuse  lag  rolling  moments.  Fuselage 
roll  angles  as  large  as  -90°  were  recorded  for  D/L  =  0.15-  At  D/L  =  0.075; 
the  fuselage  roll  angle  was  reduced  to  -30°,  and  at  a  D/L  of  0,  the 
largest  roll  angle  measured  was  -10°.  At  235  knots,  there  was  a  very  small 
change  in  roll  attitude  with  a  full  power  loss.  For  D/L=  0  and  -O.O63,  the 
fuselage  angle  of  attack  is  about  I.50.  As  seen  in  Figure  78,  rolling 
moment  is  correspondingly  small  and  independent  of  sideslip.  The  roll  dis¬ 
turbance  for  these  flight  conditions  will  be  small. 

From  Figure  20,  it  can  be  seen  that  as  D/L  increases,  the  flapping  and  lead 
angles  also  increase  following  the  loss  in  power.  Since  the  aircraft  roll 
motion  increases  with  D/L,  the  rotor  disk  lags  the  fuselage  and,  consequent¬ 
ly,  blade  flapping  also  increases. 

Effects  of  rotor  design  parameters,  D/L  =0.15 

The  effects  of  a  full  power  loss  at  D/L  =  0.15  for  all  rotor  systems  are 
shown  in  Figures  21  and  22.  The  rotor  speed  for  all  articulated  rotors 
tend.-  to  converge  to  about  70$  of  trim  speed  in  Figure  21.  The  low  Lock 
number  rotor  maintains  a  consistently  higher  rotor  speed  and  the  delta- 
three  rotor  speed  is  lowest  throughout  the  time  period.  The  severity  of 
the  disturbance  can  be  seen  in  the  yaw  rate  and  sideslip  angle  curves.  The 
peak  values  are  twice  as  large  as  for  the  part  los ^  of  power  data  shown  in 
Figure  18.  This  effect  is  consistent  with  the  twofold  increase  in  unbal- 
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anced  tail  rotor  yawing  moment  associated  with  the  full  power  loss.  Roll 
displacement  is  doubled  so  that  the  flapping  amplitude  is  markedly  in¬ 
creased  for  all  configurations.  The  fuselage  motion  of  the  delta -three 
rotor  was  so  large  that  reasonable  value  limits  of  roll  rate  (Appendix  II) 
programed  on  the  computer  were  exceeded.  The  computer  went  into  a  hold 
mode  and  the  time  history  recording  was  terminated  1.5  seconds  after  the 
loss  of  power.  This  was  also  true  for  the  zero  offset,  hingeless,  and 
floating  hub  rotors  as  seen  in  Figure  22.  Therefore,  the  severest  condi¬ 
tion  for  a  loss  of  power  is  at  D/L  =  0.15  (200  knots).  Auto-stabilization 
by  means  of  blade  pitch-motion  feedback  and  gyroscopic  damping  of  hingeless 
blades  is  not  effective  in  minimizing  the  effects  of  the  fuselage  rolling 
moment  due  to  sideslip.  Increasing  the  blade  inertia  of  the  articulated 
rotors  and  decreasing  the  rotor  propulsive  force  requirements  are  the  only 
effective  means  of  reducing  the  initial  disturbance  following  loss  of  power. 

Effects  of  rotor  design  parameters,  D/L  =  -Q.O63 

A  comparison  of  the  effects  of  a  negative  rotor  propulsive  force  for  several 
rotor  systems  is  presented  in  Figure  23.  The  responses  of  the  fully  articu¬ 
lated  rotor  (rotor  a),  delta-three  rotor  (rotor  b),  low  Lock  number  rotor 
(rotor  e),  and  hingeless  rotor  (rotor  f)are  compared  for  a  full  power  loss 
at  D/L  =  -O.O63  (235  knots).  Although  the  initial  yaw  and  the  sideslip 
disturbance  are  nearly  the  same  as  for  D/L  =  0.15  conditions,  the  roll  mo¬ 
tion  is  negligible.  As  discussed  previously,  the  fuselage  rolling  moment 
is  small  and  roughly  independent  of  sideslip  angle,  at  this  fuselage  angle 
of  attack,  so  that  all  rotor  systems  have  about  the  same  response.  Notice 
that  the  left  roll  induces  a  nose-down  pitch  motion  on  the  hingeless  rotor. 
Rotor  speed  also  decays  fastest  for  this  rotor,  and  reduces  to  70$  after  b 
seconds.  However,  the  low  Lock  number  rotor  maintains  a  rotor  speed  slight¬ 
ly  higher  than  the  standard  rotor,  and  rotor  speed  for  the  delta-three 
rotor  falls  slightly  below  the  conventional  rotor.  Since  rotor  speed  is 
the  only  quantity  that  has  changed,  the  gradual  increase  in  flapping  and 
lead-lag  amplitude  is  due  to  the  reduction  in  rotor  speed. 

Part  Power  Loss  With  Corrective  Control  Inputs,  D/L  =0.15 

The  effects  of  control  inputs  on  helicopter  and  blade  transient  response 
following  a  partial  loss  of  power  at  D/L  =  0.15  are  shown  in  Figures  2b 
through  29  for  the  following  rotor  systems:  fully  articulated  rotor  (a), 
delta-three  feedback  (b),  tip  path  plane  tilt  feedback  (c),  reduced  Lock 
number  (d),  zero  flapping  hinge  offset  (e),  and  hingeless  rotor  (f). 

It  has  been  shown  in  Figures  18  and  19  that  for  this  flight  condition  the 
critical  quantity  affecting  the  helicopter  disturbance  is  fuselage  roll  due 
to  sideslip.  Since  the  primary  cause  is  the  unbalanced  tail  rotor  thrust 
yawing  moment,  corrective  control  inputs  applied  only  to  tail  rotor  blade 
pitch  were  investigated  for  this  condition.  The  tail  rotor  control  inputs 
for  this  case  are  listed  in  Table  IV  and  consisted  of  two  ramp  inputs  of 
7.5  degrees  per  second  applied  after  a  time  delay  of  0.6  and  0.2  seconds, 

4a  and  4b  respectively,  and  a  step  input  at  0.2  second,  2B.  For  all  cases, 
the  tail  rotor  blade  pitch  angle  was  reduced  to  3  degrees. 
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The  ramp  input  at  0.6  second,  4A,  reduced  the  roll  period  and  displacement 
for  all  rotor  systems.  Ramp  and  step  inputs  at  0.2  second  reduced  roll 
angle  even  further;  and  the  step  input,  2B,  provided  the  largest  reduction 
in  roll  amplitude.  For  the  articulated  rotor  systems,  Figures  24  through 
28  show  that  reducing  the  roll  amplitude  also  decreases  the  blade  motion, 
with  a  maximum  reduction  occurring  when  the  roll  amplitude  is  minimized. 
Although  rotor  speed  for  the  articulated  rotors  is  not  significantly  af¬ 
fected  by  the  control  inputs,  the  hingeless  rotor  speed  time  history  shows 
marked  changes,  Figure  29-  The  reduction  in  aircraft  motion  caused  by  in¬ 
put  2B  decreased  the  rotor  speed-up  and  roll  motion  of  the  aircraft.  By 
reducing  the  rolling  motion  of  the  aircraft,  the  nose  down  pitch  rate 
(caused  by  gyroscopic  coupling)  was  also  decreased. 

The  strong  coupling  between  fuselage  sideslip  and  roll  is  clearly  evident 
in  these  figures.  By  reducing  tail  rotor  thrust,  the  fuselage  yawing 
moment  reverses  the  yaw  rate,  sideslip  is  rapidly  reduced  and  roll  distur¬ 
bance  is  correspondingly  reduced.  Since  the  ramp  input  with  a  0.2-second 
time  delay,  input  4b,  provided  the  least  overall  sideslip  motion,  the  ef¬ 
fects  of  the  fuselage  rolling  moment  should  be  small.  A  comparison  of 
rotor  configurations  with  this  input  will  then  permit  an  evaluation  of  the 
effects  of  rotor  parameters.  This  has  been  done  and  the  results  are  shown 
in  Figures  30  and  31* 

The  transient  response  of  the  fully  articulated  rotor  (a),  delta -three 
feedback  rotor  (b),  tip  path  plane  feedback  rotor  (c)  and  low  Lock  number 
rotor  (e)  are  compared  in  Figure  30.  A  7*5  degrees  per  second  ramp  tail 
rotor  blade  pitch  angle  input  (4b)  was  applied  0.2  second  after  the  part 
power  loss.  The  hingeless  rotor  (f)  and  zero  offset  rotor  (d)  are  compared 
with  rotor  (a)  for  the  same  control  input  in  Figure  31-  A  comparison  of 
results,  shown  in  Figure  30  indicates  that  there  is  no  significant  differ¬ 
ence  in  the  transient  response  of  these  articulated  rotors.  Rotor  speed 
varies  slightly  and  only  small  differences  occur  in  roll  displacement.  The 
low  Lock  number  rotor  rolls  slightly  less  than  the  conventional  rotor,  and 
the  introduction  of  blade  motion  feedback  produces  a  slightly  larger  roll 
angle.  Flapping  amplitude  behaves  in  the  same  way.  The  hingeless  rotor 
and  zero  offset  rotor  showed  the  only  noticeable  difference,  as  seen  in 
Figure  31*  The  zero  offset  rotor  time  history  data  was  terminated  earlier 
because  roll  velocity  exceeded  limit  values.  The  zero  offset  rotor  can, 
therefore,  be  considered  to  respond  the  most  severely  to  the  part  power 
loss  at  200  knots. 

The  hinge  less  rotor  response  is  also  relatively  severe.  Rotor  speed 
initially  falls  off  and  then  reverses  to  a  peak  value  of  about  11 6$  in 
about  3-8  seconds  after  the  power  loss.  Also,  although  the  hingeless  rotor 
rolls  slower  than  the  fully  articulated  rotor  (a),  the  peak  value  is  larger 
and  the  roll  motion  induces  a  nose  down  pitching  motion.  A  direct  compari¬ 
son  of  the  results  presented  in  Figures  30  and  31  indicates  that  blade 
inertia  and  flapping  hinge  offset  are  important  rotor  design  parameters 
which  affect  helicopter  and  rotor  blade  motion  following  loss  of  power. 

Full  Power  Loss  With  Corrective  Control  Input 

Effect  of  rotor  design  parameters  and  0.2  second  delay,  P/L  =  0.1$ 
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Although  the  transient  response  of  a  helicopter  following  a  partial  loss 
of  power  can  be  adequately  controlled  with  only  tail  rotor  corrective  con¬ 
trol  inputs,  cyclic  and  collective  controls  must  be  applied  when  a  full 
loss  of  power  occurs.  These  additional  controls  are  required  in  order  to 
recover  the  sharp  reduction  of  rotor  speed,  as  well  as  to  reduce  the  air¬ 
craft  motion. 

Since  a  full  power  loss  at  D/L  =  0.15  (200  knots)  represents  the  most 
severe  flight  condition,  several  combinations  of  control  inputs  were  inves¬ 
tigated  at  this  condition.  The  effects  of  three  control  input' patterns  on 
the  response  of  rotors  (a)  through  (f)  are  compared' in  Figures  32  through 

38. 

The  strong  effect  of  the  unbalanced  yawing  moments  following  the  full  power 
loss  is  evident  in  the  transient  roll  and  sideslip  motion  curves  for  the 
articulated  rotors  (Figures  32  through  36).  The  most  effective  means  of 
reducing  this  undesirable  roll  motion  was  found  to  be  a  step  input  on  tail 
rotor  pitch  to  0  degrees  at  0.2  second.  A  tail  rotor  step  input  at  0.6 
second  reduced  the  roll  motion  to  only  65°.  Although  roll  amplitudes  as 
large  as  100  degrees  have  been  flown  satisfactorily  by  pilots,  they  are 
considered  very  undesirable  under  unalerted  situations  such  as  this. 

Rotor  speed  was  effectively  controlled  by  reducing  collective  and  pulling 
back  on  longitudinal  cyclic  stick,  inputs  5B,  12B,  and  12C.  Input  5B  is 
not  the  optimum  control  input  for  this  flight  condition  and  would  require 
further  corrective  action  by  the  pilot  within  2  seconds  (Figures  32  through 
37)-  Input  patterns  12B  and  12C  were  found  to  be  the  most  promising 
(Figure  37)-  Figure  37  also  shows  the  input  of  rotor  control  (collective 
and  cyclic  pitch)  can  be  delayed  up  to  four  seconds,  following  a  loss  in 
power,  and  still  control  rotor  speed  satisfactorily.  Since  all  articulated 
rotor  configurations  respond  in  a  similar  manner  to  control  inputs  IB,  5B, 
and  6b,  it  is  expected  they  would  all  respond  similarly  to  input  12B  and 
12C.  The  traces  with  input  12B  terminated  at  1.4  seconds  because  qw2  ex¬ 
ceeded  its  built-in  limit  (Appendix  II) .  The  value  of  qw2  also  became  very 
close  to  its  built-in  limit  for  control  input  12C.  This  is  not  a  physical 
limit  but  is  an  expected  maximum  limit  based  on  test  data.  Further  rescaling 
would  have  prevented  qW2  from  becoming  too  large. 

Improvement  in  roll  response  was  achieved  by  applying  right  lateral  stick, 
input  oB  (Figures  32  through  36).  The  total  effect  of  this  input  is  an 
order  of  magnitude  reduction  in  the  roll  motion  of  an  articulated  rotor 
following  a  full  power  loss  at  D/L  =  0.15.  In  all  cases,  vibratory  flap¬ 
ping  amplitude  was  initially  reduced  by  this  input,  but  was  then  increased 
prior  to  the  termination  of  the  run. 

The  most  effective  initial  control  input  for  the  hingeless  rotor,  input  8b, 
is  almost  identical  to  input  6b  for  the  articulated  rotors.  The  only  dif¬ 
ference  is  that  collective  pitch  was  reduced  to  8°  instead  of  2°.  With 
this  input,  the  aircraft  rolled  left  only  to  about  -5°  and  then  recovered 
(Figure  38).  Rotor  speed  also  recovered  after  a  mild  droop  to  about  92$. 
The  aircraft,  however,  pitched  nose  down  to  -10°.  This  nose  down  tendency 
contrasts  with  the  sharp  nose  up  effect  caused  by  input  9B.  For  this  input 
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collective  was  reduced  only  to  12°  and  longitudinal  control  was  changed  by 

3-75°. 

The  transient  response  associated  with  control  input  6b  for  the  articulated 
rotors  and  8b  for  the  hingeless  rotor  has  been  redrawn  in  Figures  39  and 
40  in  order  to  obtain  a  relative  comparison  for  all  systems.  Lock  number 
and  flapping  hinge  offset  consistently  stand  out  as  important  design  param¬ 
eters.  The  delta -three  rotor  shows  a  smaller  recovery  of  rotor  speed  than 
the  other  articulated  rotor  systems  for  the  same  control  input.  This  is 
due  to  the  washout  of  control  input  with  delta-three;  therefore,  it  would 
require  more  control  input  to  achieve  the  same  amount  of  rotor  control  as 
the  other  articulated  rotors.  Other  than  the  direct  effect  of  the  control 
input  on  blade  flapping,  there  is  little  or  no  relative  difference,  with 
or  without  control  inputs,  in  the  response  of  the  articulates  rotor  systems. 

Effect  of  0.6  second  time  delay,  fully  articulated  rotor,  P/L;  0 

The  effect  of  applying  control  inputs  following  a  full  power  loss  with  a 
0.6-second  time  delay  was  investigated  for  the  fully  articulated  rotor  (a) 
at  D/L  =  0  (2Q0  knots),  and  the  results  are  shown  in  Figure  4l.  Although 
the  tail  rotor  input,  3A,  ultimately  reduced  sideslip  to  zero,  the  initial 
effects  of  rolling  moment  produced  a  left  roll.  As  sideslip  was  reduced  to 
zero  no  fuselage  moment  was  available  to  restore  the  aircraft  to  a  level 
attitude.  However,  the  aircraft  maintained  a  reasonably  level  attitude 
and  some  rotor  speed  was  recovered  with  the  full  control  input,  7A.  Since 
this  input  represents  a  typical  pilot  applied  input,  the  pilot  should  ex¬ 
perience  no  difficulty  in  recovering  from  a  full  loss  of  power  at  this 
flight  condition,  D/L  =  0. 

Effect  of  rotor  design  parameters  at  239  knots,  D/L  =  -0.063 

The  effect  of  control  input  time  delay  has  also  been  investigated  at  235 
knots,  and  the  results  are  presented  in  Figures  42  through  47.  The  effects 
of  ramp  tail  rotor  control  inputs  with  a  0.2-second  time  delay,  input  3B, 
and  a  0.6-second  time  delay,  input  3A,  are  compared,  along  with  a  full  con¬ 
trol  input,  7B,  in  Figures  42  through  4-7  for  four  rotor  systems.  These 
rotor  systems  include  the  fully  articulated  rotor  (a)  in  Figure  42,  the 
delta-three  rotor  (b)  in  Figure  43,  the  reduced  Lock  number  rotor  (e)  in 
Figure  44,  and  the  hingeless  rotor  (f)  in  Figure  45.  A  relative  comparison 
of  the  effects  of  control  input  7B  is  presented  in  Figure  4-6  for  all  of 
these  rotors. 

Since  the  roll  disturbance  is  small  following  a  full  power  loss  at  the 
flight  condition  D/L  =  -O.063,  tail  rotor  control  inputs  were  applied  to 
reduce  the  yaw  rate  and  sideslip  motion.  The  rapid  application  of  tail 
rotor  corrective  controls,  0.2-second  time  delay,  effectively  reduced  both 
the  yaw  and  sideslip  disturbance  and  the  flapping  amplitude.  As  the  tail 
rotor  thrust  is  reduced,  a  left  rolling  moment  is  generated.  This  happens 
because  the  main  rotor  rolling  moment  is  no  longer  balanced,  so  the 
aircraft  rolls  left.  If  the  tail  rotor  control  input  is  delayed  0.6  sec¬ 
ond,  the  roll  motion  is  smaller.  This  is  true  for  all  configurat  ons.  The 
additional  lateral  stick  motion  in  input  7E  reduced  this  roll  motion,  the 
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yaw  rate,  sideslip  angle,  and  the  rotor  speed  decay  for  all  the  articulated 
rotors . 

The  gyroscopic  coupling  of  the  hingeless  rotor,  however,  caused  a  nose 
down  pitching  motion  following  the  tail  rotor  control  inputs  as  shown  in 
Figure  4 5.  Although  the  roll  response  was  slightly  reduced  when  the  delay 
time  was  increased  to  0.6  second,  the  pitching  motion  was  not  significantly 
affected.  The  lateral  stick  motion  of  input  10B  reduced  sideslip.  Yaw  rate 
was  reduced  only  at  first  but  then  was  increased  to  the  right.  The  pitch 
response  was  aggravated  by  this  control  input. 

The  nose  down  motion  of  the  hingeless  rotor  is  contrasted  with  the  slight 
pitch  up  tendency  of  the  fully  articulated  rotor  systems  in  Figure  46.  The 
hingeless  rotor  pitched  nose  down  to  about  30°  in  2.5  seconds  and  because 
of  this  would  require  further  pilot  corrective  action  to  maintain  attitude 
control.  The  delta-three  rotor  systems  response  was  milder  than  the  other 
articulated  rotor  systems.  This  was  due  to  the  washout  of  control  inputs 
caused  by  the  delta -three  coupling. 

For  the  articulated  rotor  at  a  D/L  =  -O.063,  the  most  promising  control 
input  was  found  to  be  a  reduction  in  main  rotor  collective  pitch  to  2°  and 
longitudinal  change  in  cyclic  pitch  of  -5-5°>  input  13B  (Figure  47).  Since 
the  low  Lock  number  and  articulated  rotor  responded  in  a  similar  manner  to 
control  input  3A,  3B,  and  10B,  the  low  Lock  number  rotor  would  be  expected 
to  respond  similarly  to  input  13B.  For  this  high  speed  flight  condition, 
although  roll  motion  control  presents  no  problem,  application  of  cor¬ 
rective  tail  rotor  control  inputs  to  recover  the  yaw  and  sideslip  motion 
must  be  accompanied  by  right  lateral  stick  motion.  Aft  stick  and  reduced 
collective  will  recover  the  rotor  speed  droop.  The  disturbance  following 
a  full  power  loss  at  235  knots,  D/L  =  -O.O63,  is  mild,  and  the  pilot 
should  have  no  difficulty  in  controlling  the  aircraft.  A  full  power  loss 
at  D/L  =  0.15,  however,  produces  a  severe  disturbance  and  an  automatic 
tail  rotor  control  system  would  probably  be  required  in  order  to  control 
the  aircraft  fuselage  motion. 

C.  EFFECT  OF  ROTOR  DESIGN  PARAMETERS,  ROTOR 
PROPULSIVE  FORCE,  AND  CORRECTIVE  CONTROL 
INPUTS  ON  BLADE  BENDING  MOMENTS 


Time  history  data  of  blade  bending  moments  were  recorded  for  all  rotor 
systems  except  for  the  floating  hub  rotor.  Data  were  obtained  for  27  com¬ 
binations  of  rotor  configuration  and  control  inputs.  Since  the  most  severe 
aircraft  and  blade  disturbance  occurred  at  the  D/L  =  0.15  condition,  thir¬ 
teen  combinations  are  included  at  this  flight  condition.  Included  is  each 
of  the  six  configurations  with  no  control  inputs,  and  data  for  three  con¬ 
trol  inputs  apiece  for  the  articulated  rotor  (a)  and  the  delta-three  rotor 
(b).  Typical  time  history  data  recorded  for  this  study  are  presented  in 
Figure  48.  Multiple  moment  data  are  shown  for  the  fully  articulatea  rotor 
(a)  following  a  full  power  loss  at  200  knots,  (D/L  s  0.15).  It  is  seen  in 
this  figure  that  flapwise,  chordwise,  and  torsional  moment  data  were  all 
recorded  at  three  radial  stations.  Stations  between  the  nodal  points  of 
the  first  three  flapwise  bending  modes  were  selected  for  this  study.  The 
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radial  stations  for  the  articulated  rotor  are  4l,  56,  and  71  percent  of  the 
rotor  radius.  The  hingeless  rotor  was  investigated  at  22. 5,  45,  and  70 
percent.  In  addition,  chordwise  bending  moment  data  in  the  vicinity  of  the 
lag  damper,  r/R  =  12$,  were  recorded  for  the  articulated  rotors.  The  re¬ 
corded  moments  can  be  converted  to  stress  by  means  of  the  section  modulus 
data  listed  in  Table  VI.  Since  the  neutral  axis  is  not  equidistant  from 
the  leading  and  trailing  edges,  two  values  are  given  for  the  chordwise  sec¬ 
tion  modulus. 

Effect  of  Rotor  Design  Parameters  and  Drag-to-Lift  Ratio  on  Trim  Bending 
Moments 


The  trim  bending  moments  for  all  the  study  conditions  listed  in  Table  II 
have  been  recorded  and  the  results  are  presented  in  Figures  49  through  52. 
The  one-half  peak-to-peak  vibratory  moment  is  plotted  as  a  function  of  the 
corresponding  steady  moment  to  constitute  a  bending  moment  diagram  for  each 
of  the  four  drag-to-lift  ratios. 

The  data  for  flapwise  bending  of  the  articulated  rotors  at  r/R=5 6$  (  Figure 
49)  indicate  that  the  low  Lock  number  rotor  consistently  has  the  lowest 
steady  flapwise  moment,  and  has  a  vibratory  flapwise  moment  slightly  higher 
than  the  other  systems  for  all  drag-to-lift  ratios.  For  the  same  forward 
speed,  200  knots,  the  vibratory  and  steady  moments  generally  decrease  as 
the  rotor  propulsive  force  requirements  are  reduced  from  D/L  =  0.15  to  0. 
This  reduction,  however,  is  due  in  part  to  a  reduction  of  blade  twist  which 
is  -4°  at  D/L  =0.15  and  0.075  and  0°  at  D/L  =0  and  -O.O63.  The  differ¬ 
ence  in  the  vibratory  bending  moment  among  rotor  systems,  however,  is  small. 
At  D/L  =  0.15,  the  difference  between  the  vibratory  flapwise  moment  for  the 
low  Lock  number,  highest  value,  and  the  delta-three  rotor,  lowest  value, 
is  equivalent  to  only  ±  600  psi  of  blade  stress. 

The  chordwise  steady  bending  moments  for  all  rotor  systems  are  negative 
except  for  the  zero  offset  rotor,  as  seen  in  Figure  50.  The  zero  offset 
rotor  also  has  the  smallest  vibratory  chordwise  bending  moment;  since  the 
chordwise  moments  presented  in  this  figure  correspond  to  less  than  2000  psi 
steady  and  vibratory  stress,  the  values  are  not  critical  for  the  trim  con¬ 
ditions  shown. 

The  torsional  vibratory  moment  is  also  smaller  at  D/L  =  0  than  at  D/L  = 

0.15,  as  seen  in  Figure  51*  For  all  drag-to-lift  ratios  the  steady  and 
vibratory  torsional  moment  is  the  smallest  for  the  fully  articulated  rotor, 
and  generally,  the  delta-three  rotor  experiences  the  largest  vibratory 

torsional  moment. 

The  hingeless  rotor  trim  bending  moments  are  shown  in  Figure  52.  An  in¬ 
board  reference  station  at  r/R  s  22-5 $  was  selected  since  the  maximum 
moments  occur  in  the  inboard  section  of  a  hingeless  rotor.  The  negative 
steady  chordwise  bending  moments  could  have  been  minimized  if  the  rotor 
blades  had  been  swept  aft. 
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Effect  of  Rotor  Design  Parameters,  Rotor  Propulsive 
Force  and  Control  Inputs  on  Blade  Vibratory  and 
Peak  Moment  Time  History  for  Full  Power  Loss 

Time  history  of  the  peak  and  vibratory  blade  bending  moment  envelopes  are 
presented  in  Figures  53  through  58  for  each  of  the  three  spanwise  stations. 
These  moment  envelopes  will  provide  a  thorough  appreciation  of  the  effects 
of  rotor  parameters  and  control  inputs  on  blade  bending  moments. 

The  peak  moment  is  the  variation  of  the  maximum  negative  or  positive  moment 
recorded  as  a  function  of  time.  Therefore,  the  envelope  of  the  peak  moment 
is  defined  by  these  curves.  Also  shown  is  the  one-half  peak-to-peak  moment 
envelope,  called  the  vibratory  moment. 

Effect  of  rotor  design  parameters,  P/L  ■  0.15,  QQ  control  input 

Various  articulated  configurations  are  compared  in  Figure  53  at  a  rotor 
lift-to-drag  ratio  of  0.15  for  a  full  power  loss.  All  configurations  shown 
have  about  the  same  flapwise  and  chordwise  moment  envelope  following  a  full 
power  loss,  except  in  the  region  from  1  to  3  seconds,  where  there  is  an  in¬ 
crease  in  flapwise,  chordwise,  and  torsional  steady  and  vibratory  moments 
with  configurations  (a)  and  (c).  This  moment  increase  is  due  to  retreating 
blade  stall,  which  is  shown  in  Figure  48  as  peaks  in  torsional  moment  at 
1.35  and  1.75  seconds.  The  termination  of  configuration  (b)  at  1.7  sec¬ 
onds  was  due  to  fuselage  roll  rate  exceeding  its  built  in  value  of  1  rad./ 
sec.  Retreating  blade  stall  causes  an  increase  in  flapwise  moment  on  the 
outboard  section  of  the  blade,  while  the  maximum  chordwise  moment  is  in¬ 
creased  along  the  entire  blade. 

The  maximum  flapwise  and  chordwise  vibratory  moments  recorded  for  config¬ 
urations  (a),  (b),  and  (c)  were  converted  to  stress,  using  Table  V,  and  cor¬ 
respond  to  vibratory  stresses  below  the  5>500  psi  allowable  vibratory 
stress  of  aluminum  for  effective  infinite  life.  Therefore,  this  maneuver 
does  not  result  in  a  vibratory  or  peak  bending  moment  which  would  cause 
damage  to  the  blades.  Blade  stress  is  used  here  as  a  criteria  to  determine 
the  severity  of  the  maneuver.  The  moment  loading  in  this  study  can  be 
used  for  analyzing  other  rotor  systems,  provided  the  El  ratio  and  blade 
bending  natural  frequencies  of  the  other  system  is  the  same  as  that  used 
in  this  study.  In  addition,  the  stress  produced  by  these  moments  will  be 
the  same  if  the  i/C  values  of  the  two  systems  are  the  same. 

The  zero  offset  rotor  system  and  low  Lock  number  rotors  are  compared  to 
articulated  rotor  in  Figure  54  and  at  a  rotor  drag-to-lift  ratio  of  0.15 
for  a  full  loss  in  driving  torque.  Reducing  Lock  number  increases  the 
chordwise  and  flapwise  peak  and  vibratory  moments  on  the  blade  following  a 
power  failure.  The  flapwise  vibratory  moment  at  the  outboard  radial  sta¬ 
tion,  where  the  maximum  moment  occurs,  corresponds  to  a  stress  level  of 
about  10,000  psi.  At  this  stress  level  the  blade  could  withstand  more  than 
3,000,000  cycles,  therefore,  it  is  not  considered  a  serious  condition  for 
this  blade  design.  The  most  pronounced  effect  of  the  zero  offset  rotor, d, 
is  the  reduced  trim  chordwise  peak  and  vibratory  moments.  Following  the 
power  failure,  the  flapwise  and  chordwise  moments  build  up  in  the  same  manner 
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as  the  articulated  rotor.  The  termination  of  the  zero  offset  rotor  after 
one  second  was  due  to  the  fuselage  roll  rate  exceeding  its  built  in  limit. 

It  is  of  interest  to  note  even  with  such  a  -violent  fuselage  motion  the  flap- 
wise  and  chordwise  blade  moments  were  not  exceedingly  large;  near  the  end 
of  the  trace  the  torsional  moment  did  show  an  increase.  A  typical  time 
history  variation  of  the  chordwise  moment  in  the  vicinity  of  the  lag  damper 
(12$  radial  station)  is  shown  in  Figure  48  for  the  articulated  rotor. 
Throughout  this  study,  this  moment  remained  small. 

Effect  of  Control  inputs,  fully  articulated  rotor,  P/L  =  0.1$ 

The  effect  of  control  inputs  on  blade  bending  moments  is  shown  in  Figure  55 
for  the  articulated  rotor  at  a  rotor  D/L  of  0.15  for  a  full  loss  in  power. 

By  reducing  fuselage  yaw  and  roll  with  a  step  input  on  tail  rotor  collective 
pitch,  the  maximum  torsional  moment  was  reduced.  The  flapwise  and  chordwise 
vibratory  moments  recorded  correspond  to  stresses  below  the  allowable  5>500 
psi  for  infinite  life  and,  therefore,  are  not  serious.  Even  the  chordwise 
and  flapwise  bending  moments  produced  by  control  input  5B  and  6b  correspond 
to  stresses  below  10,000  psi  and  are  not  serious. 

The  initial  effect  of  the  main  rotor  control  input  is  to  reduce  chordwise 
and  torsional  bending  moments.  This  is  due  to  unloading  the  blades  with 
the  reduction  in  collective  pitch.  After  one  second  all  bending  moments 
increase  as  the  rotor  tilts  back  and  rotor  flapping  increases. 

Effect  of  Control  inputs,  fully  articulated  rotor,  D/L  =-0.063 

The  effect  of  control  inputs  on  the  articulated  rotor  blade  moments  at  a 
rotor  drag-to-lift  ratio  of-O.063  is  shown  in  Figure  5 6.  Introduction  of 
tail  rotor  pitch,  to  reduce  fuselage  sideslip,  decreases  the  flapwise  and 
chordwise  bending  moments,  but  has  very  little  effect  on  the  torsional 
moment.  Rotor  inputs  11B  initially  decrease  the  torsional  and  fiapwise 
vibratory  moments.  At  0.5  to  1.0  seconds,  the  vibratory  and  peak  moments 
increase  with  rotor  inputs.  As  in  the  previous  case,  these  moments  will 
not  damage  the  blades. 

Effect  of  rotor  propulsive  force 

The  hingeless  rotor  bending  moments  did  not  show  an  appreciable  change  fol¬ 
lowing  a  full  power  loss  until  fuselage  motion  built  up  (Figures  57  and  58). 
The  termination  of  both  traces  prior  to  5  seconds  was  due  to  the  fuselage 
roll  rate  exceeding  1  rad. /sec.  The  sudden  build  up  of  blade  bending 
moments  in  Figure  58  was  due  to  a  sudden  increase  in  fuselage  roll  rate 
caused  by  the  rigid  rotor  gyroscopic  moment.  As  the  main  rotor  slows  down, 
the  aircraft  pitches  nose  down.  This,  in  turr,  causes  a  violent  negative 
gyroscopic  rolling  moment.  If  automatic  stabilization  had  been  used,  this 
build  up  would  probably  have  been  prevented. 

The  moments  presented  here  are  those  for  an  aluminum  blade.  The  moments 
would  increase  if  steel  were  used.unless  the  El  distribution  and  blade  bend¬ 
ing  natural  frequencies  were  kept  the  same  as  that  of  the  rotor  blade  used 
in  this  study. 
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TABLE  IV 

CONTROL  INPUT  PATTERNS 


Ramp  Inputs  At  7*5°  Per  Second 


Input 

Number* 


Lower 

Limit 


Lower  Limit 


Input  Increment 


*  A.  Input  Applied  0.6  Second  After  Power  Loss. 

B.  Input  Applied  0.2  Second  After  Power  Loss. 

C.  Main  Rotor  Input  Applied  4.0  Seconds  After  Power  Loss 
And  Tail  Rotor  Input  0.2  Second  After  Power  Loss. 


STEP  INPUT  PATTERN 
| Trim  Value 


Control 

Input 


Lower  Limit 


Control 

Input 


Time 

Delay  " 

RAMP  INPUT  PATTERN 

Trim  Value 


7 • 5°/ sec , 


Time 

Delay' 


Time 


Lower  Limit 
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TABLE  V 


TRIM  CONDITIONS 


Configuration 

D/L 

dt 

0C 

Als 

®ls 

af 

mm 

Fully  Articulated 

0.15 

7.4 

16.5 

-4.9 

11.0 

-4 

1.8 

Rotor  (a) 

0.075 

6.1 

14.3 

-4.0 

10.5 

-1.3 

1.2 

0.0 

4.9 

6.9 

-3-0 

6.7 

1.5 

1.0 

-0.0b3 

4.0 

4 .9- 

-3-2 

3-5 

1-5 

1.0 

Delta -3  Feedback 

0.15 

7.8 

17.8 

-6.5 

10.3 

-4.5 

1.0  ' 

Rotor  (b) 

0.0 

4.1 

8.1 

-4.4 

6.4 

1.5 

.0 

-O.O63 

4.9 

7-7 

-3-7 

4.5 

2.0 

1.0 

Tip  Path  Plane 

0.15 

7.5 

lo .  6 

-0.2 

11.7 

-4.2 

1.0 

Tilt  Feedback 

0.0 

4.3 

7-2 

-4.0 

7.2 

1.6 

1.0 

Rotor  (c) 

Zero  Offset 

0.15 

9.0 

17.6 

-7.075 

10.2 

-6.6 

1.2 

Rotor  (d) 

0.0 

5.0 

8.4 

-6.8 

8.9 

2.0 

1.0 

Low  Lock  Number 

0.15 

7.9 

16.7 

-3.6 

11.7 

-4.5 

1.0 

Rotor  (e) 

0.075 

6.4 

14.5 

-2.9 

11.2 

-1.5 

1.0 

0.0 

4.4 

7.0 

-2.4 

6.8 

1.5 

1.0 

-O.O63 

4.5 

5-3 

-2.4 

4.8 

2.0 

1.0 

Hingeless 

0.15 

7-4 

15.6 

-2.6 

10.0 

-5.0 

1.0 

Rotor  (f) 

0.075 

7-1 

12.8 

-3.2 

8.1 

-2.0 

1.0 

0.0 

6.5 

c  .8 

-5.9 

5.6 

•5 

1.0 

-O.O63 

0.1 

3-8 

-0.9 

2  9' 

1.5 

1.0 

Floating  Hub 

0.15 

8.85 

14.6 

-3-9 

6.5 

-6.5 

1.0 

Rotor  (g) 

0.075 

6.0 

12.7 

-4.1 

8.4 

-1.0 

1.0 

0.0 

3-7 

5.6 

-2.2 

4.4 

2-5 
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TABLE  VI 


SPAR  SECTION  MODULUS 


ARTICULATED  ROTOR 

SYSTEMS 

Radius  Rotor,  r/R 

Flapwise 

Chordwise 

Torsional 
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Two  Seconds  After  Power  Loss 
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Figure  6.  Effect  of  Lock  Number  on  Flapping 
Disturbance  Following  Power  Loss 
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Figure  7.  Effect  of  Lock  Number  on  Blade  Lead-Lag 
Motion  Following  Power  Loss 
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Figure  9«  Effect  of  Pitch-Blade  Motion  Feedback  on 
Rotor  Speed  Two  Seconds  After  Power  Loss 
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Figure  10.  Effect  of  Pitch-Blade  Motion  Feedback  on 
Flapping  Disturbance  Following  Power  Loss 
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Figure  11.  Effect  of  Pitch-Blade  Motion  Feedback  on 

Blade  Lead-Lag  Motion  Following  Power  Loss 
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Figure  l4.  Effect  of  Flapping  Hinge  Offset  on  Flapping 
Disturbance  Following  Power  Loss 
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Figure  15 .  Effect  of  Flapping  Hinge  Offset  on  Blade 
Lead-lag  Motion  Following  Power  Loss 


42 


BLADE  FLAP  ANGLE,  /9.DEG. 
LEAD  ANGLE,  8,  DEG. 


10 


ROTOR  a 


TIME,  SEC. 


Figure  l8.  Concluded 


b6 


PERCENT  ROTOR 


BLADE  FLAP  ANGLE,  fi,  DEG. 
LEAD  ANGLE  8,  DEG. 


mi 

81 

IBI 

IBI 

■Mil 

IB 

Bl 

WM 

Bl 

IBI 

ii 

flAW] 

p 

Jill 

Ml 

IBBP 

IIVII 

IVIi 

m 

IM1 

mi 

nil 

SIB! 

Bl  Bl 

rJM 

TAv 

Bl 

mmu 

Bl 

t  PI 

VMM 

Bl 

jSa 

Inn 

ill 

|| 

IB 

B 

m 

|| 

MB 

■MM 

IHI 

imwM 

IaIPJU 

BMI 

IBI 

Ifll 

IB 

III 

IB 

UVS 

Ml 

IB 

HWM 

Ml 

Bl 

Mil 

m 

■ 

IBB 

iasai 

IMS 

Ml 

Uwa 

m 

Ml 

IB 

vwml 

Ml 

BBto 

Brivk’ 

Ml 

MrM 

SB 

kwi 

B 

IM 

S3 

BE 

g  MM 

in 

■■IB 

1 

n 

IbI 

IB 

Ml: 

Ml 

Ml 

mmE 

■Bl 

■M 

IB 

ip 

Bri 

BMM 

IBB 

Inn 

B 

$ 

B 

1 

B 

IB 

IBB 

b 

mam 

MM 

p 

II 

— 1 — 

-4 - H 

■pi 

■Ml 

ip 

In 

IB 

■Ml 

IBB 

Mil 

Rbbi 

IHPf 

■ 
■  ■■ 

B 

■ 

1MI 

IBI 

■1  II 

IB 

■B" | 

IB 

II 

llll 

Immi 

■  ■■ 

Ml 

111! 

■ 

B 

B 

■1 

IBI 

II 

Ell 

| 

MM 

M1W 

3S 

S 

Ml 

M 

B 

Ml 

■Ml 

IBI 

IB 

|KS2! 

llll 

ji 

■ia 

Ml 

M 

B 

B 

II  1 

Ml 

Ml 

■■1 

ibi 

IBI 

I! 

Bl 

IBI 

0  1.0  2.0  3.0  4.0 


TIME,  SEC. 

Figure  19 .  Concluded 

48 


FUSELAGE  SIDESLIP  PITCH  RATE  YAW  RATE,  t 


20 


0 

•20 


TIME,  SEC. 


Figure  20.  Helicopter  and  Roto±  Blade  Motion  Following  Full 
Power  Loss,  No  Corrective  Control  Inputs,  Fully 
Articulated  Rotor  (e/R  =  0.03*0 
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PERCENT  ROTOR 


Figure  21.  Helicopter  and  Rotor  Blade  Motion  Following  Full 

Power  Loss,  D/L  =  0.15,  No  Corrective  Control  Inputs 
Various  Rotor  Configurations 
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Figure  21.  Concluded 


Figure  22.  Helicopter  and  Rotor  Blade  Motion  Following  Full 

Power  Loss,  D/L  a  0.15,  No  Corrective  Control  Inputs, 
Various  Rotor  Configurations 
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Helicopter  and  Rotor  Blade  Motion  Following  Full 
Power  Loss,  D/L  =  -  O.O63,  No  Corrective  Control 
Inputs,  Various  Rotor  Configurations 
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Figure  2k.  Helicopter  and  Rotor  Blade  Motion  Following  Part 
Power  Loss,  D/L  z  0.15,  With  Corrective  Control 
Inputs,  Fully  Articulated  Rotor  (e/R  =  0.034) 
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Figure  25.  Helicopter  and  Rotor  Blade  Motion  Following  Part 
Power  Loss,  D/L  =  0.15,  With  Corrective  Control 
Inputs,  Rotor  With  Blade  Pitch-Flap  Feedback 
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Figure  26.  Concluded 
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Figure  27.  Helicopter  and  Rotor  Blade  Motion  Following  Part 

Power  Loss,  D/L  =  0.15,  With  Corrective  Control  Inputs, 
Rotor  With  Reduced  Lock  Number 
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Figure  28.  Concluded 
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Figure  29.  Helicopter  and  Rotor  Blade  Motion  Following  Part 

Power  Loss,  D/L  =  0.15?  With  Corrective  Control  Inputs, 
Hingeless  Rotor 
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Figure  31*  Helicopter  and  Rotor  Blade  Motion  Following  Part 
Power  Loss,  D/L  =  0.15,  With  Corrective  Control 
Input  4b,  Various  Rotor  Configurations 
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Figure  32.  Helicopter  and  Rotor  Blade  Motion  Following  Full 

Power  Loss,  D/L  z  0.15,  With  Corrective  Control  Inputs 
Fully  Articulated  Rotor  (e/R  =  0.03*0 
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Figure  33.  Helicopter  and  Rotor  Blade  Motion  Following  Full 

Power  Loss,  D/L  r  0.15,  With  Corrective  Control  Input 
Rotor  With  Blade  Pitch-Flap  Feedback 
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Figure  34.  Concluded 
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Figure  35-  Helicopter  and  Rotor  Blade  Motion  Following  Full 

Power  Loss,  D/L  =  0.15,  With  Corrective  Control  Inputs 
Rotor  With  Reduced  Lock  Number 
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Figure  36.  Concluded 
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Figure  37.  Helicopter  and  Rotor  Blade  Motion  Following  Full  Power  Loss, 
D/L  =  0.15,  With  Corrective  Control  Inputs,  Fully  Articulated 
Rotor  (e/R  =  0.034) 
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Figure  39-  Helicopter  and  Rotor  Blade  Motion  Following  Full 

Power  Loss,  D/L  =  0.15,  With  Corrective  Control  Input  — 

oB,  Various  Rotor  Configurations 
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Figure  40.  Helicopter  and  Rotor  Blade  Motion  Following  Full  -v 

Power  Loss,  D/L  =  0.15,  With  Corrective  Control  Input 
6b,  Various  Rotor  Configurations 
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Figure  4l.  Concluded 
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Figure  42.  Helicopter  and  Rotor  Blade  Motion  Following  Full 

Power  Loss,  D/L  =  -  O.O63,  With  Corrective  Control 
Inputs,  Fully  Articulated  Rotor  (e/R  =  0.034) 
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Figure  44.  Helicopter  and  Rotor  Blade  Motion  Following  Full 

Power  Loss,  D/L  =  -  O.O63,  With  Corrective  Control 
Inputs,  Rotor  With  Reduced  Lock  Number 
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Helicopter  and  Rotor  Blade  Motion  Following  Full 
Power  Loss,  D/L  =  -  O.O63,  With  Corrective  Control 
Inputs,  Hinge less  Rotor 
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Figure  46.  Helicopter  and  Rotor  Blade  Motion  Following  Full 

Power  Loss,  D/L  r  -  O.O63,  With  Corrective  Control 
Input  7B,  Various  Rotor  Configurations 
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Figure  kj .  Helicopter  and  Rotor  Blade  Motion  Following  Full 

Power  Loss,  D/L  =  -  O.O63,  With  Corrective  Control 
Inputs,  Fully  Articulated  Rotor  (e/R  =  0.034) 
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Figure  48.  Fully  Articulated  Rotor  (e/R  =  0.034),  Full  Power 
Loss,  D/L  z  0.15,  No  Corrective  Control  Input, 
Multiple  Bending  Moment  Data 
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Figure  48.  Continued 
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Figure  48.  Concluded 
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Figure  50.  Trim  Chordwise  Moment  Characteristics  (r/R  =  56$), 

Various  Rotor  Configurations  and  Drag-To-Lift  Ratios 
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Figure  52.  Hingeless  Rotor  Trim  Bending  Moment  Characteristics 
(r/R  =  22.5$),  Various  Drag-To-Lift  Ratios 
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Figure  55-  Rotor  Blade  Vibratory  and  Peak  Moment  Variation, 

Full  Power  Loss,  l)/L  s  0.15,  Various  Corrective  Con¬ 
trol  Inputs,  Fully  Articulated  Rotor  (e/R  =  0.034) 
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CONCLUSIONS 


1.  A  full  loss  of  main  rotor  power  can  seriously  affect  the  safety  of  a 
200  -  235  knot  helicopter,  unless  there  is  a  rapid  and  appropriate  applica¬ 
tion  of  tail  rotor  control.  For  the  S-61F  helicopter,  used  in  this  study, 
operating  at  a  rotor  propulsive  force -to-lift  ratio,  D/L  =  0.15,  tail 
rotor  corrective  control  action  was  necessary  in  less  than  0.2  second  to 
retain  attitude  control.  The  loss  of  power  is  less  critical  with  reduced 
rotor  propulsive  force.  With  all  propulsion  provided  by  auxiliary  propul¬ 
sion  devices,  corrective  control  input  delay  times  of  0.6  second  for  both 
main  rotor  and  tail  rotor  control  inputs  can  be  tolerated.  Rotor  speed 
recovery  was  achieved  using  collective  and  cyclic  pitch  corrective  control 
with  time  delays  up  to  4.0  seconds  (Figure  59)* 

2.  High  speed  helicopters  should  be  configured  with  large  vertical  and 
horizontal  tail  surfaces  to  insure  a  stable  platform  in  the  event  of  an 
abrupt  loss  of  power.  The  aircraft  response  is  du?  primarily  to  the  speci¬ 
fic  fuselage  aerodynamic  characteristics  associated  with  the  sideslip  which 
is  created  by  the  unbalanced  yawing  moment.  The  weathercock  stability 
provided  by  a  large  vertical  tail  will  limit  the  sideslip  angle  and  hence 
reduce  the  resulting  rolling  moment  due  to  the  fuselage  and  rotor  dihedral 
effect  (Figure  60).  The  beneficial  effect  of  a  large  horizontal  tail  is 
clearly  evident  in  the  small  fuselage  pitch  disturbance  and  stable  motion 
observed  for  the  S-61F  base  helicopter  throughout  this  study. 

3.  Flapping  hinge  offset  and  roJ  or  blade  inertia  (Lock  number)  are  impor¬ 
tant  design  parameters  affecting  helicopter  and  blade  motion  and  loading 
following  a  loss  of  power.  A  rotor  with  a  moderate  flapping  hinge  offset 
ratio,  3  to  5  percent,  experienced  both  a  smaller  roll  disturbance  than  the 
zero  offset  rotor  and  a  smaller  blade  moment  excursion  than  the  larger  ef¬ 
fective  offset  rotor  (hingeless).  A  reduction  in  Lock  number  results  in 
reduced  fuselage  rolling  and  blade  flapping  disturbance.  However,  for  the 
conditions  investigated, the  higher  inertia  of  the  rotor  blades  increases 
the  tendency  for  the  blade  to  stall  due  to  roll  rate,  and  results  in  higher 
values  of  blade  bending  moments  (Figure  6l). 

4.  Blade  pitch-flap  feedback  (Delta-three  r  30°)  and  blade  pitch-tip  path 
plane  feedback  are  not  effective  concepts  for  reducing  the  helicopter  roll 
motion  following  a  loss  of  power.  Rotor  blade  stall  tendency  (caused  by 
the  larger  roll  rate)  and  the  corresponding  blade  bending  moments  were  in¬ 
creased  by  Delta-three  and  tip  path  plane  tilt  feedback.  Below  stall,  both 
concepts  increased  this  torsional  vibratory  and  peak  bending  moment 
(Figure  62) . 

5.  The  gyroscopic  effects  of  the  hingeless  rotor  system  restrain  the  ini¬ 
tial  roll  motion  following  power  loss,  but,  without  stability  augmentation, 
the  pitch-roll  coupling  aggravates  the  long  term  roll  and  pitch  motion 
(Figure  63) . 

6.  Blade  tip-tail  cone  clearance  will  not  be  a  problem  for  rotor  fuselage 
arrangements  typical  of  the  configuration  used  in  this  study.  The  maximum 
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down  flapping  angle  referenced  to  the  shaft  over  the  tail  cone  would  provide 
more  than  adequate  margin  for  all  reasonable  recovery  maneuvers  (Figure  64). 

7*  The  loss  of  altitude  becomes  more  pronounced  as  the  rotor  propulsive 
force  increases  following  a  loss  in  power.  This  loss  in  altitude  is  due  in 
a  large  part  to  pitch  and  roll  attitude  build  up  (Figure  65). 

8.  The  marked  increase  in  bending  moments  generated  by  the  motion  follow¬ 
ing  the  power  loss  can  be  diminished  by  reducing  the  resulting  roll  rate 
with  corrective  tail  rotor  inputs  and  also  by  unloading  the  blades.  The 
bending  moments  increase  as  the  rotor  angle  of  attack  increases  with  cyclic 
pitch  control  inputs.  The  values  of  the  bending  moments  recorded  during 
this  study,  under  fixed  control  inputs,  did  not  produce  blade  stresses 
which  would  result  in  a  loss  of  the  aircraft  for  the  duration  of  time  his¬ 
tories  shown  (Figure  66).  The  bending  moments  for  the  hingeless  rotor 
tend  to  be  higher  than  for  the  other  rotors  and  should  be  considered  in  the 
design  of  the  aircraft. 
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Figure  62.  Effect  of  Blade  Pitch -Flap  Coupling  on  Helicopter 
Motion,  Rotor  Speed,  and  Torsional  Bending  Moment 
Following  Full  Power  Loss  at  200  Knots 
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Figure  63.  Response  of  the  Hingeless  and  Articulated  Rotor 
Following  a  Part  Power  Loss  at  200  Knots 


133 


BLADE  DEFLECTION  BLADE  FLAPPING  ANGLE 

»  IN.  0 lif-o,  r/R>i  o)>  DEG, 


ALTITUDE  CHANGE,  AZ,  FT. 


Figure  65.  Change  in  Altitude  Following  Loss  of  Main  Rotor 
Power,  Fully  Articulated  Rotor  (e/R  =  0.03^-) 
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Figure  66.  Effect  of  Corrective  Control  Inputs  With  0.2 

Second  Time  Delay  on  Rotor  Blade  Bending  Moments 
Following  Full  Loss  of  Power  at  200  Knots,  Fully 
Articulated  Rotor  (e/R  =  0.03^) 


136 


RECOMMENDATIONS 


1.  This  exploratory  study  has  disclosed  several  significant  trends  con¬ 
cerning  the  effects  of  operating  conditions,  airframe  aerodynamics,  and 
corrective  control  inputs  on  helicopter  and  blade  motion  following  a  loss 
of  power  at  high  speeds,  all  of  which  should  be  explored  in  more  detail. 

The  fuselage  aerodynamic  characteristics,  in  particular,  were  found  to  be 
significant  due  to  the  dynamic  pressures  associated  with  the  high  airspeeds. 
For  the  flight  conditions  investigated,  the  rotor  sustained  all  the  air¬ 
craft  gross  weight  and  was  operating  near  the  stall  limit.  Consideration 
should  be  given  to  alleviating  rotor  stall  by  means  of  wings  and  flaps.  In 
addition,  auxiliary  stabilization  devices,  pilot  reaction  characteristics, 
extended  time  delays,  and  coupling  of  rotor  and  auxiliary  propulsion  have  a 
profound  effect  on  the  helicopter  and  blade  motion  and  on  the  corrective 
control  inputs  following  a  loss  of  power.  Since  these  considerations  were 
beyond  the  scope  of  this  study,  it  is  recommended  that  this  study  be  con¬ 
tinued  to  examine  in  greater  detail  the  effects  of  the  following  areas: 

.  wing  and  flap  arrangements 

.  fuselage  aerodynamic  and  stability  characteristics 
.  auxiliary  stabilization  equipment 
.  control  inputs  from  a  simulated  pilot  including 
extended  time  delays 

.  coupled  power  requirements  of  rotor  and  auxiliary 
propulsion  systems 

2.  This  study  has  re-emphasized  the  importance  of  the  time  delay  before 
corrective  control  inputs  are  applied.  The  problem  facing  the  helicopter 
pilot  experiencing  a  loss  of  engine  power  is  complex,  both  because  the  cues 
are  complex  (noise,  accelerations,  vibrations,  angular  rates)  and  because 
the  decisions  he  must  make  involve  not  only  "what"  but  "how  much".  Experi¬ 
mental  data  on  pilot  reactions  to  an  unalerted  power  loss  in  flight  are 
almost  nonexistent.  Research  on  the  three  following  parts  of  the  pilot  de¬ 
lay  time  problem  is  therefore  recommended: 

.  tests  to  establish  the  cues  that  can  be  detected 
by  a  pilot 

.  a  specific  study  of  helicopter  pilot  reaction  times 
.  flight  tests  with  simulated  emergencies 

3.  Although  the  method  developed  in  this  study  has  shown  good  correlation 
with  flight  test  data  at  143  knots,  further  correlation  of  fuselage  motion, 
blade  motion,  and  loading  are  required  at  higher  airspeeds.  Flight  test 
data  for  maneuvers  at  high  airspeeds  will  be  obtained  in  the  S-6lF  (NH-3A) 
Maneuvering  Flight  Test  Program  to  be  conducted  under  Army/Navy  sponsorship 
in  1966.  The  data  obtained  from  this  program  should  be  correlated  with  the 
analytical  methods  developed  in  this  present  study. 
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APPENDIX  I 


ROTOR  SYSTEM  AND  AIRCRAFT  PHYSICAL 
PROPERTIES 

ROTOR  SYSTEM  PHYSICAL  PROPERTIES 


The  physical  properties  of  all  of  the  rotor  systems  are  listed  in  Table  VII, 
and  the  spanwise  physical  properties  and  bending  mode  shapes  of  their  blades 
are  presented  in  Figures  67  through  72.  The  aerodynamic  characteristics  of 
the  blade  airfoil  sections  are  presented  in  the  Development  of  Computer 
Program  Section.  The  S-61R  (CH-3C)  rotor  was  used  in  the  correlation 
phase,  and  rotors  (a)  through  (g)  were  used  in  the  analytical  phase  of  this 
study.  The  S-6lR  main  rotor  blade,  with  appropriate  modifications,  is  com¬ 
mon  to  all  of  the  articulated  rotor  systems,  rotors  (a)  through  (e).  A 
scaled-down  version  of  the  blades  used  on  the  C^-R  compound  helicopter  of 
Reference  6  was  used  for  the  rotor  systems,  with  blades  rigidly  attached  to 
the  hub,  rotor  systems  (f)  and  (g).  These  rotors  were  sized  to  have  the 
same  loading  as  the  articulated  rotors. 

Since  the  blade  torsional  moment  of  inertia  is  determined  experimentally 
for  the  whole  blade  including  the  cuff,  the  spanwise  distribution  is  not 
known.  For  this  study,  the  torsional  mode  shapes  were  calculated  by  assum¬ 
ing  a  constant  spanwise  distribution  of  torsional  moment  of  inertia,  which 
was  derived  from  the  experimental  value  for  the  whole  blade. 

The  bending  mode  shapes,  their  derivatives,  and  the  spanwise  physical  prop¬ 
erties  of  the  CH-3C  rotor  blade  are  presented  in  Figure  67.  The  da^a  in¬ 
clude  three  flap-wise  mode  shapes  with  their  first  and  second  derivatives, 
two  chordwise  mode  shapes  with  their  first  and  second  derivatives,  and  one 
torsional  mode  shape  with  its  first  derivative.  The  spanwise  physical 
properties  include  flaj wise,  chordwise,  and  torsion  moments  of  inertia  as 
well  as  section  area  and  mass  distribution. 

The  blades  used  in  rotor  systems  (a),  (b),  and  (c)  are  identical,  except 
for  twist,  to  the  S-6lR  blades.  Since  the  tip  speed  is  633  fps  compared 
with  660  fps  for  the  S-6lR,  the  mode  shapes  will  have  slight  differences. 
The  mode  shapes  and  spanwise  physical  properties  for  these  rotors  are 
presented  in  Figure  b8. 

The  zero  offset  rotor  blade,  rotor  system  (d),  was  also  derived  from  the 
S-61R.  The  root  structure  and  mass  were  extended  to  the  center  of  rota¬ 
tion.  The  resulting  structural  and  physical  characteristics  are  shown  in 
Figure  69. 

A  reduced  Lock  number  for  rotor  system  (e)  was  obtained  by  adding  a  non- 
structural  tip  weight  of  2b.  5  pounds  to  the  blades  of  rotor  system  (a). 

With  this  weight,  a  26$  reduction  in  Lock  number  from  10.62  to  7-88,  was 
achieved.  The  spanwise  data  for  these  blades  a~e  presented  in  Figure  70. 

The  torsional  moment  of  inertia  was  increased  to  account  for  the  blade  ti] 
mass . 


ibO 


The  torsional  moment  of  inertia  for  the  hingeless  blades,  shown  in  Figures 
71  and  72,  was  empirically  derived  from  production  blade  characteristics 
and  was  assumed  constant  along  the  span.  Since  the  end  conditions  are  dif¬ 
ferent,  the  mode  shapes  for  the  fixed  hub  blades,  shown  in  Figure  71>  dif¬ 
fer  from  the  mode  shapes  for  the  floating  hub  blades  presented  in  Figure 
72.  The  selection  of  the  mode  shapes  is  discussed  in  the  Development  of 
Computer  Program  Section  of  this  report.  In  both  cases,  the  blade  was 
cantilevered  from  an  infinitely  rigid  hub.  The  hub  extended  outboard  to  a 
9 (j>  spanwise  station. 

The  torsional  mode  shape  for  all  rotors  includes  the  effect  of  control  sys¬ 
tem  flexibility.  The  root  spring  constant  for  the  S-6lR  push  rod-swash- 
plate  system  was  determined  experimentally.  The  root  spring  constant  for 
the  hingeless  blades,  however,  was  increased  to  maintain  the  same  root  to 
blade  spring  ratio  as  the  S-61R. 

AIRCRAFT  PHYSICAL  PROPERTIES 


The  physical  properties  of  the  S-61R  (CH-3C)  aircraft  (Figure  73)>  and  the 
S-6lF  (NH-3A)  high  speed  research  aircraft  (Figure  jb),  used  in  this  study 
are  presented  in  Table  VIII.  The  loading  conditions  and  mass  moments  of 
inertia  of  the  S-6lF  aircraft,  including  the  effective  inertia  of  the  rotor 
system,  were  kept  constant  for  all  rotor  systems. 

The  weight  and  balance  data  are  based  on  a  low  loading  condition.  Although 
the  center  of  gravity  position  is  well  aft,  the  large  stabilizing  surfaces 
make  the  aircraft  stable  for  this  condition. 

The  aerodynamic  characteristics  of  both  aircraft,  obtained  from  wind  tunnel 
model  tests,  are  presented  in  Figures  75  and.  76.  This  data  was  transformed 
from  the  wind  axis  system  to  the  body  axis  system  for  use  in  this  program 
and  resulted  in  the  curves  shown  in  Figures  77  and  78.  Since  all  lift 
forces  were  assumed  to  be  generated  entirely  by  the  rotor  system,  the  fuse¬ 
lage  and  stabilizer  lifting  forces  and  pitching  moments  were  adjusted  for 
each  rotor  system  and  speed  to  meet  the  rotor  operating  conditions  listed 
in  Table  III. 

A  conventional  helicopter  control  system  of  longitudinal  and  lateral  cyclic, 
collective,  and  tail  rotor  blade  pitch  was  assumed.  The  controls  were  not 
coupled  with  each  other,  nor  with  any  auxiliary  surfaces,  such  as  elevators 
or  rudders.  The  elevator  was  assumed  to  be  trim  adjustable  to  provide 
adequate  cyclic  stick  control  range. 


TABLE  VII 


ROTOR  PHYSICAL  PROPERTIES 


Rotor  System 

CH-3C 

ra 

d 

e 

t,s 

Radius,  ft. 

31.0 

31.0 

31.0 

31.0 

28.7 

Offset 

Ratio,  5 

.0339 

.0339 

0 

.0339 

.0900 

Chord,  in. 

18.25 

18.25 

18.25 

18.25 

29.52 

b 

5 

6 

6 

6 

4 

*1 

-8° 

0°, -4° 

0°,  -4° 

0°,  -4° 

0°, -4° 

{i0R,  ft. /sec. 

660 

633 

633 

633 

633 

WB,  lb. 

275.5 

275.5 

332.5 

300.0 

399-0 

MB(rcg)'  slug-ft. 

92.5 

92.5  • 

103.0 

114.1 

110.2 

p 

IB,  slug-ft. 

1800 

1800 

2002 

2430 

1918 

y 

10.62 

10.62 

9.56 

7.88 

11.86 

Ihub>  slug-ft.2 

573-3 

573-3 

568.8 

573-3 

625.9 

Ar 

.07661 

.07661 

.07661 

.07661 

.07661 

( e/ r ) LD 

.776 

.776 

1.000 

.776 

-- 

e  ♦  e' 

•  2339 

•  2339 

•  2339 

•  2339 

•  2339 

w  w^/ St  o 

2.678 

2.689 

2.568 

2.843 

1.300 

0 

Ci 

3 

4.889 

4.943 

4.507 

5.303 

3.830 

w  w^/ St  o 

7.631 

7.772 

6.959 

8.305 

8.270 

<•)  Vj/Jl  o 

3.454 

3.530 

3.310 

2.213 

1.433 

C 

rd 

o 

9.105 

9.390 

8.490 

8.440 

7.518 

0 

CJ 

Q> 

3 

8.070 

8.397 

8.397 

8.057 

8.154 
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TABLE  VII  (Continued) 


Rotor  System 

CH-3C 

a,b,  c 

d 

e 

ffg 

Kr  X  10 "6 

1.200 

1.200 

1.200 

1.200 

5-300 

Co/lB 

.632 

.632 

.632 

.632 

-- 

AR  Blade 

20.  4 

20.4 

20.4 

20.4 

11.7 

*Spar  Cavity 

Area,  in. 

11.40 

11.40 

11.40 

11.40 

38.49 

Ip,  lbs. -ft. - 
sec .  “ 

.672 

.672 

.672 

•  938 

2.932 

ca 


M 


TABLE  VIII 

AIRCRAFT  PHYSICAL  PROPERTIES 


CH-3C 

S-61F 

*  GW,  lb. 

16,000 

11,900 

**  lx,  slugs  -  ft/ 

8,650 

10,180 

**  Iy,  slugs  -  ft. 

52,500 

53,100 

*#  Iz,  slugs  -  ft.2 

47,300 

51,200 

C.G.  J 

Water  Line,  in. 

153.6 

164.0 

Fuselage  Station,  in. 

280.0 

278.0 

Main  Rotor 

Water  Line,  in. 

233.0 

233.0 

Fuselage  Station,  in. 

267.5 

267.5 

Tail  Rotor 

Water  Line,  in. 

225.0 

233.5 

Fuselage  Station,  in. 

709.5 

705.0 

Aerodynamic  Data  Reference 

Water  Line,  in. 

127.4 

144.0 

Fuselage  Station,  in. 

275.0 

275.0 

f  Horizontal  Tail 

Water  Line,  in. 

255.0 

ISS.O 

Fuselage  Station,  in. 

709.5 

701.50 

Area,  ft.^ 

27.0 

76.0 

Vertical  Fin  I 

Water  Line,  in. 

200.0 

186.8 

Fuselage  Station,  in. 

672.0 

686.2 

Area,  ft.2 

24.0 

42.5 

Shaft  Incidence,  deg.,  Forward 

4.0 

3-0 

*  Including  Rotor  System. 

**  Including  the  Effective 

Inertia 

of  Rotor  System. 
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FLAPWISE  NORMAUZED  BLADE 


NORMALIZED  BLADE  RADIUS,  */R 

Figure  67.  CH-3C  Rotor  Blade  Spanwise  Structural  and 
Physical  Properties 
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BLADE  SECTION  FLAPWISE  MOMENT  OF  INERTIA,  l„.  IN.4  SECOND  DERIVATIVE  OF  FLAPWISE 

NORMALIZED  BLADE  DEFLECTION.  T, 


Figure  67.  Continued 
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Figure  67 .  Continued 
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NORMALIZED  BLADE  RADIUS,  x/R 


Figure  68.  Articulated  Rotor  Blade  Spanwise  Structural 
and  Physical  Properties 
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NORMALIZED  BLADE  RADIUS,  x/R 


.1  .2  3  .4  .5  .6  .7 

NORMALIZED  BLADE  RADIUS,  x/R 


Figure  68.  Continued 
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VISE  NORMALIZED  BLADE 
DEFLECTION,  X, 
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NORMALIZED  BLADE  RADIUS,  x/R  * 

Figure  68.  Continued 
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NORMALIZED  BLADE  RADIUS,  x/R 


Figure  68.  Continued 


Figure  68.  Concluded 
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NORMALIZED  BLADE  RADIUS,  x/R 


Figure  69.  Zero  Offset  Rotor  Blade  Spanwise  Structural 
and  Physical  Properties 
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Figure  69 
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D  BLADE  RADIUS,  x/R 


Continued. 


CHORDWISE  NORMALIZED  BLADE 


Figure  69.  Continued 
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NORMALIZED  BLADE  RADIUS,  x/H 


Figure  69.  Continued 
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Figure  69.  Continued 
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NORMALIZED  BLADE  RADIUS,  x/R 


Figure  70.  Reduced  Lock  Number  Rotor  Blade  Spanwise 
Structural  and  Physical  Properties 
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Figure  70.  Continued 
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NORMALIZED  BLADE  RADIUS,  x/R 


Figure  70-  Continued. 


FIRST  DERIVATIVE  OF  NORMALIZED  NORMALIZED  BLADE  TORSIONAL 

BLADE  TORSIONAL  DEFLECTION,  YJ  DEFLECTION, 


Figure  70.  Concluded 
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NORMALIZED  BLADE  RADIUS,  x/R 

Figure  71*  Hingeless  Blade  Spanwise  Structural  and 
Physical  Properties 


SECTION  CHORDWtSE  MOMENT  OF  INERTIA,  1M,  IN.4  SECOND  DERIVATIVE  OF 
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FIRST  DERIVATIVE  OF  NORMALIZED  NORMALIZED  BLADE  TORSIONAL 


Figure  71.  Continued 
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NORMALIZED  BLADE  RADIUS,  x/R 


Figure  71*  Concluded 
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Fi/yire  7  .  Floating  Hub  Rotor  Blade  Spanvise  Structural 
and  Physical  Properties 
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SECOND  DERIVATIVE  OF  FLAPWISE 

BLADE  SECTION  FLAPWISE  MOMENT  OF  INERTIA,  I„.  IN.4  NORMALIZED  BLADE  DEFLECTION.  X 
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FIRST  DERIVATIVE  OF  CHORDWISE  CHORDWISE  NORMALIZED  BLADE 

NORMALIZED  BLADE  DEFLECTION,  Yj  DEFLECTION,  X, 
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Figure  72.  Continued 
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Figure  72.  Continued 
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Figure  72.  Concluded 
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Figure  7^.  S-6lF  (NH-3A)  Aircraft 


PITCHING  MOMENT 


Figure  75.  S-61R  (CH-3C)  Wind  Tunnel  Data  (Wind  Axis  System) 
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Figure  75*  Concluded 
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Figure  76.  Concluded 
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Figure  78*  S-61F  (NH-3A)  Wind  Tunnel  Data  (Body  Axis  System) 
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Figure  Jd.  Concluded 
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EQUATIONS  OF  MOTION  AND  DESCRIPTION  OF 
FLEXIBLE  BLADE  AND  HYBRID  COMPUTER 
SIMULATION 


EQUATIONS  OF  MOTION 


In  anticipation  of  large  transient  aircraft  motion  following  an  abrupt  loss 
of  power,  equations  of  motion  of  the  fuselage  were  developed  with  no  small 
angle  assumptions.  The  rotor  forces  and  moments  were  initially  referenced 
to  a  shaft  coordinate  system  and  the  fuselage  wind  tunnel  data  were  initial¬ 
ly  referenced  to  a  "wind  tunnel"  axis  system.  All  forces  and  moments  were 
then  transfered  to  the  body  axis  system  for  a  final  solution  of  the  general 
equations  of  motion. 

The  origin  of  the  body  axis  system  is  located  at  the  aircraft  center  of 
gravity.  The  X-axis  is  positive  forward  along  the  fuselage  horizontal  ref¬ 
erence  line  in  the  plane  of  symmetry.  The  Z-axis  is  also  in  the  plane  of 
symmetry  but  perpendicular  to  the  X-axis  and  positive  down.  The  Y-axis  is 
perpendicular  to  the  X-Z  plane  and  positive  to  the  right. 

The  inertial  loading  of  the  system  is  treated  in  the  conventional  manner. 
However,  since  the  principal  axes  are  very  close  to  the  aircraft  reference 
axes,  the  products  of  inertia  are  considered  reasonably  small  and,  there¬ 
fore,  are  not  considered.  The  classic  Euler  equations  were  used  to  calcu¬ 
late  motion  in  the  space  axis  system  for  the  inertia  loads.  The  effect  of 
rotor  downwash  on  fuselage  and  horizontal  tail  angle  of  attack  is  included 
in  this  analysis. 

The  general  body  axis  equations  of  motion  are  as  follows: 

Along  Body  X  axis: 

m  (^Xf  *  c*uZf”  ruYf)  =  sin  ^  +  Pr°PuXsive  Force 

Along  Body  Y  axis: 

m  (uYf-+  rux  -  p  ugf)  s  Yj®  ♦  Yff  -f  Tt  ♦  mg  cos  9  sin  $ 

Along  Body  Z  axis: 

m  (uZf  +  p  Uyf-  quxf)  s  Zmr  -  Lff  ♦  mg  cos  9  cos  ^ 

About  the  Body  X  axis: 

p  Ix  -  qr  (Iy-Iz)  8  hMRxMR  ♦  4lR  +  hTR  ^T  +  **ff  +  ha  Yff 
About  the  Body  Y  axis: 

4  Iy  -  rP  (XZ  -  XX^  s  XMR  ZMR  -  hMR  XMR  + 
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Mj^jr  ♦  la  Lff  ♦  haDff+  Mdamp  ♦  Mdownwasn 


About  the  Body  Z  axis: 

r  Iz  -  pq  (lx  "  Iy)  =  "  ^MR  ^MR  ♦  OMR  “ 

Ijr  Tp  Ndamp 


The  transformations  of  the  various  applied  forces  and  angular  relationships 
can  be  found  in  Appendix  III  along  with  further  explanation  and  details  of 
the  above  equations.  A  more  simple  Bailey  analysis  was  used  in  the  calcu¬ 
lation  of  tail  rotor  thrust. 

DESCRIPTION  OF  FLEXIBLE  BLADE  SIMULATION 


The  vibratory  motion  of  the  flexible  rotor  blades  was  simulated  by  a  super¬ 
position  of  normal  modes.  This  was  accomplished  by  first  deriving  the  dif¬ 
ferential  equations  of  motion  for  the  lateral  and  torsional  deformation  of 
a  linearly  twisted,  mass  balanced  helicopter  rotor  blade  operating  in  a 
maneuvering  flight  condition.  The  equations  were  initially  derived  from  a 
consideration  of  the  equilibrium  of  aerodynamic,  dynamic,  and  elastic 
moments  acting  on  a  given  blade  section.  The  resulting  moment  equations 
were  then  differentiated  to  be  expressed  in  terms  of  the  local  blade  load¬ 
ings.  These  differential  equations  were  then  solved  by  simulating  the 
blades  with  a  series  of  normal  modes. 

The  flapwise,  chordwise,  and  torsional  elastic  deflections  were  expressed 
as  finite  series  summations  of  assumed  radial  shape  functions  (mode  shapes), 
with  each  suitably  scaled  by  time -dependent  generalized  coordinates.  The 
number  of  terms  retained  constitutes  three  flapwise,  two  chordwise,  and  one 
torsional  mode  shape.  The  mode  shapes  used  in  this  study  correspond  to  the 
natural  vibratory  mode  shapes  of  an  untwisted  blade  mounted  on  a  fixed  hub 
and  operating  at  zero  pitch,  flap,  and  lag  angles  and  at  the  trim  rotor 
speed. 

The  substitution  of  these  mode  shapes  permits  the  three  equations  of  flap- 
wise,  chordwise,  and  torsional  loading  to  be  expanded  into  a  system  of  dif¬ 
ferential  equations,  with  the  generalized  coordinates  as  unknowns.  Further, 
because  of  the  orthogonality  properties  of  the  particular  mode  shapes 
selected,  modal  dynamic  coupling  terms  have  been  eliminated.  Additional 
simplifications  were  achieved  by  replacing  the  structural  and  centrifugal 
stiffness  terms  with  natural  frequency  terms.  Chordwise  structural  damping 
was  also  included  in  the  analysis.  An  equivalence  of  5$  structural  damping 
was  simulated  by  2.5$  viscous  damping,  Reference  4. 

Rotor  equations  of  motion  were  developed  and  programed  so  that  only  minor 
modifications  were  required  for  each  configuration.  For  the  hingeless 
rotor,  several  program  expressions  were  altered;  two  degrees  of  freedom, 
flapping  and  lead-lag  were  removed;  and  flapping  angle  was  replaced  by  a 
fixed  coning  angle  of  4°.  Since  the  blades  are  rigidly  attached  to  the  hub, 
the  rotor  speed  degree  of  freedom  equation  was  simplified  by  using  only 
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aerodynamic  torque.  This  procedure  also  eliminated,  any  algebraic  instabil¬ 
ity  from  dynamic  coupling. 

The  floating  hub  rotor  was  simulated  by  a  four  blade  rotor  with  two  ortho¬ 
gonal  teetering  rotors.  Initially,  all  elastic  deflections  were  measured 
relative  to  a  plane  which  was  tangent  to  the  hub.  Due  to  an  undesirable 
mathematic  instability,  the  method  was  revised  so  that  the  average  tilt  of 
the  rotor  disk  was  used  as  a  plane  of  reference.  This  procedure  required 
the  introduction  of  pinned-free  as  well  as  fixed-free  bending  modes,  as 
shown  in  Figure  72.  Due  to  the  complexity  of  this  change,  the  floating  hub 
analysis  was  restricted  to  inelastic  degrees  of  freedom.  Thus,  only  heli¬ 
copter  and  blade  motion  data  were  obtained  for  this  configuration. 

Some  of  the  assumptions  on  which  the  derivation  of  the  flexible  blade  equa¬ 
tions  of  motion  was  based  are  as  follows: 

1.  The  blade  section  center  of  gravity  is  coincident  with  the  intersection 
of  the  blade  elastic  axis  and  the  principal  axis  of  the  section. 

2.  All  elastic  deflections  and  rotations  are  small. 

3.  Blade  flapping  and  lagging  angles  are  small,  so  that  sin  (  )  = 

(  )  and  cos  (  )=  l-(  f  2. 

4.  Fuselage  angular  and  linear  velocities  are  small  compared  with  rotor 
angular  velocity  and  tip  speed,  respectively. 

5.  Quasi -steady  aerodynamic  theory  is  applicable,  and  rotor  inflow  is 
uniform  across  the  disk. 

6.  Radial  flow  effects  are  negligible. 

7.  Apparent  mass  aerodynamic  forces  are  negligible. 

8.  Blade  twist  is  both  small  and  linear  with  radius. 

9-  The  flapping  and  lagging  hinges  are  offset  from  the  center  of  rotation 
by  less  than  10$  of  the  rotor  radius. 

10.  The  blade  flapping  and  lagging  hinges  are  coincident  for  the  fully 
articulated  rotors. 

11.  The  change  in  radial  position  of  a  blade  mass  element  due  to  elastic 
deformations  is  negligible  when  compared  to  rigid  body  blade  motions 
for  the  articulated  and  floating  hub  rotors. 

12.  The  thickness -to-chord  and  the  chord-to-radius  ratios  of  the  blade 
are  small. 

13-  Pitch-flap  coupling  affects  only  the  blade  aerodynamics. 
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14.  The  lag  hinge  offset  for  the  zero  offset  rotor  can  he  represented  by 
an  equivalent  linear  spring. 

These  simplifying  assumptions  are  standard  for  rotor  blade  aeroelastic 
analyses  and  are  considered  adequate  to  qualitatively  describe  the  effects 
of  a  power  loss  at  high  speeds.  Correlation  studies  for  the  Sikorsky  blade 
aeroelastic  analysis  have  been  discussed  in  References  2,  and  6.  Com¬ 
parisons  of  calculated  data  using  constant  inflow  with  flight  test  data  in¬ 
dicated  that  calculated  performance  and  one -half  peak-to-peak  flapwise  blade 
stresses  gave  good  agreement  with  flight -measured  values  at  l40  knots. 
However,  as  concluded  in  Reference  2,  a  better  definition  of  induced  veloc¬ 
ity  was  required  to  predict  higher  harmonic  air  loads  and  blade  stresses 
with  greater  accuracy  at  all  airspeeds.  The  Cornell  Aeronautical  Labora¬ 
tory  (CAL)  variable  inflow  analysis,  as  used  in  Reference  1,  did  not  pro¬ 
duce  peak-to-peak  stress  values  matching  the  good  correlation  obtained 
using  a  constant  inflow  analysis  at  the  high  speed  conditions.  For  high 
speed  flight,  then,  the  assumption  of  constant  inflow  was  considered  valid. 

DESCRIPTION  OF  HYBRID  COMPUTER  SIMULATION 


A  hybrid  computer,  which  included  both  a  digital  and  analog  computer,  was 
used  to  solve  the  coupled  equations  of  motion,  in  order  to  take  full  ad¬ 
vantage  of  the  more  desirable  features  of  each  system.  The  hybrid  facil¬ 
ity  located  at  the  United  Aircraft  Research  Laboratories  was  used  for  this 
study.  This  facility  consists  of  a  general  purpose  digital  computer, 
Digital  Equipment  Corporation  PDP-1,  and  a  general  purpose  analog  computer, 
Beckman  Ease  2133.  The  two  computers  are  logically  linked  for  mode  control 
and  arithmetically  linked  by  20  high  speed  conversion  channels  in  each 
direction. 

The  PDP-1  computer,  used  to  solve  the  helicopter  equations  of  motion  and 
rotor  dynamics,  is  a  single  address,  single  instruction,  stored  program 
machine  with  a  word  length  of  18  bits.  The  memory  unit  contains  a  capacity 
for  storing  1^,336  words  with  a  memory  cycle  time  of  5  microseconds.  The 
computer  instructions  that  are  available  permit  a  wide  range  of  logic  opera¬ 
tions  and  fixed  point  arithmetic.  The  2133  analog  computer  was  used  to 
provide  automatic  controls  for  commanding  trim,  pilot  inputs  after  power 
loss,  and  output  data  for  time  history  recording.  Some  of  the  computer 
equipment  includes  100  operational  amplifiers  and  150  coefficient  potenti¬ 
ometers.  In  addition  to  the  linear  equipment,  the  computer  also  contains 
automatic  potentiometer  setting  and  component  read-out  equipment. 

The  dynamic  response  of  the  helicopter  and  the  rotor  blades  can  be  separ¬ 
ated  into  two  frequency  bands.  The  six  elastic,  flapping  and  lead-lag 
degrees  of  freedom  constitute  a  high  frequency  range.  The  rotor  speed  and 
six  helicopter  degrees  of  freedom  appear  in  the  lower  frequency  range.  The 
high  frequency  degrees  of  freedom  were  integrated  on  the  digital  computer 
by  means  of  a  Z  transform  technique,  Reference  3*  Also,  a  computer  time 
scale  of  l/lOO  of  real  time  was  selected  to  insure  conversion  of  the  higher 
bending  modes.  Simple  Euler  integration  was  found  adequate,  in  this  time 
frame,  for  the  lower  frequency  degrees  of  freedom. 
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The  rotor  and  helicopter  equations  of  motion  are  solved  by  the  digital  com¬ 
puter  as  a  three  loop  system.  A  simplified  flow  diagram  of  the  computer 
program  is  presented  in  Figure  79.  The  program  includes  aerodynamic  calcu¬ 
lations  for  each  of  ten  radial  stations  along  the  blade,  forces,  moments, 
and  dynamics  for  up  to  six  blades,  and  forces,  moments,  and  dynamics  for 
rotor  speed  and  helicopter  motion.  Although  the  program  is  in  closed  loop 
form,  it  starts  with  the  calculation  of  forces  and  moments  and  concludes 
with  calculation  of  accelerations.  In  order  to  synchronize  ohe  digital 
computer  with  the  analog  operations,  a  digital  clock  maintains  the  program 
cycle  at  a  one  second  time  interval. 

At  the  start  of  each  clock  interval,  the  program  converts  data  from  the 
analog  computer  into  the  digital  computer.  The  data  converted  includes 
pilot  control  inputs  and  helicopter  attitudes  and  velocities.  During  the 
computer  time  allocated  for  data  conversion,  the  computer  shares  data  trans¬ 
fer  operations  with  the  calculation  of  fuselage  aerodynamics.  The  body 
axis  component  of  velocity  is  used  to  calculate  forward  speed,  angle  of 
attack,  and  sideslip  angle.  With  the  aid  of  the  fuselage  wind  tunnel  data 
(Figures  77  and  78),  the  fuselage  forces  and  moments  are  determined  and  then 
the  tail  rotor  thrust  is  calculated.  Square  root  operations  in  these  cal¬ 
culations  are  accomplished  by  four  iterations  of  the  Newton  Rhapson  square 
root  technique,  and  all  required  arctangents  are  computed  utilizing  a  three 
term  Chebyshev  polynomial. 

Preliminary  data  are  then  calculated  for  use  in  the  individual  blade  loop. 
These  include  flapping  and  lead-lag  integration  coefficients  and  a  term 
for  elastic  torsion  acceleration.  Also,  the  gravitational  force  vector, 
weight,  is  resolved  into  the  fuselage  and  shaft  axes  along  with  the  air¬ 
craft  translation  and  angular  velocities  and  accelerations.  Following  the 
calculation  of  main  rotor  inflow,  the  rotor  angular  velocity  is  calculated 
with  provisions  for  forcing  the  engine -supplied  torque  to  zero  or  one  half 
value  for  full  and  part  power  loss.  A  sudden  loss  of  power  was  simulated 
by  an  exponential  time  decay  with  a  time  constant  of  0.25  second. 

The  loop,  in  which  the  calculations  for  each  blade  are  performed,  is 
started  by  the  calculation  of  azimuth  angle  corresponding  to  the  current 
blade  position.  The  local  blade  angle  associated  with  specified  control 
inputs  is  then  determined,  along  with  the  blade  segment  velocities,  angle 
of  attack,  shears,  and  moments. 

Composite  0012  airfoil  data  in  the  form  of  lift,  drag,  and  pitching  moment 
coefficient  as  a  function  of  local  angle  of  attack  and  Mach  number,  shown 
in  Figure  80,  are  used  to  calculate  the  aerodynamic  forces  and  moments  for 
each  blade  segment.  The  calculations  for  each  blade  segment  are  repeated 
until  the  forces  and  moments  are  evaluated  at  the  center  of  ten  equally 
spaced  segments,  starting  at  blade  radial  station  0.23  and  extending  to 
the  tip. 

When  the  data  for  one  blade  have  been  determined,  the  blade  flapping,  lag¬ 
ging  and  elastic  accelerations,  velocities,  and  displacements  are  then  cal¬ 
culated.  The  accelerations  are  integrated  by  means  of  difference  equations 
which  were  derived  from  a  Z  transform  equivalence  of  a  second  order 
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undamped,  system.  Reference  3* 

With  the  fade  flapping  and  displacements  known,  the  deflection  of  the 
blade  tip,  measured  relative  to  a  plane  normal  to  the  shaft,  is  then  calcu¬ 
lated.  In  this  way  the  blade  tip  clearance  over  the  tail  cone  can  be  es¬ 
tablished.  The  tip  deflection  is  computed  for  any  blade  within  5°  azimuth 
location  of  the  tail.  The  deflection  value  in  the  time  history  traces  is 
changed  only  when  the  blade  enters  the  angular  section  of  interest.  The 
individual  blade  calculations  are  completed  with  the  calculation  of  the 
blade  root  shear  forces  and  bending  moments.  When  the  calculation  of  the 
data  for  all  the  blades  has  been  completed,  the  total  force  and  moments  are 
calculated. 

With  the  rotor  calculations  completed,  the  forces  and  moments  generated  by 
the  fuselage,  main  rotor,  tail  rotor,  and  auxiliary  jet  engines  are  summed 
in  the  general  equations  of  motion  to  determine  the  aircraft  accelerations, 
velocities,  and  displacements.  The  program  then  returns  to  the  outer  loop 
beginning  the  calculation  sequence  again  to  complete  the  clock  interval. 

If  the  aircraft  is  in  a  trimming  mode,  the  controls  are  adjusted  to  reduce 
the  aircraft  accelerations  to  zero  and  to  bring  the  forward  speed  to  its 
trim  value.  At  the  same  time  pertinent  data  are  transmitted  to  the  six- 
channel  Sanborn  recorders  for  time  history  trace. 

Due  to  the  fact  that  the  PDP-1  computer  operates  with  fixed  point  arith¬ 
metic,  and  the  digital-analog  converters  require  fixed  point  data  of  a 
specified  scale,  predetermined  scale  factors  had  to  be  introduced  into  the 
computer  program.  Scaling  factors  were  selected  to  permit  as  large  a  range 
as  possible  in  parameters  without  severely  restricting  accuracy.  These 
scale  factors,  Table  IX,  were  based  on  flight  test  information  of  maximum 
values  expected  to  be  incurred  on  present  day  helicopters  and  do  not  neces¬ 
sarily  represent  physical  limits.  These  scaling  factors,  therefore,  lim¬ 
ited  the  range  of  some  of  the  parameters.  The  computer  is  programed  to  go 
into  a  hold  mode  when  these  limits  are  exceeded.  Minor  scale  changes  were 
made  during  this  study  to  insure  that  sufficient  history  data  were  acquired 
for  analysis. 


196 


TABLE  IX 


SCALE  FACTORS 


Rotor  Systems 

CH-3C,a,b,c, e 

d 

f 

g 

9  deg. 

±180 

±l8o 

±180 

±l8o 

9  f,  deg. /sec. 

±1 

±  1 

±  1 

±1 

^  Xf'  ^Yf'^Zf 

±1 

±  1 

±  1 

±1 

Axf’AvAZf  1/sec- 

±1 

±  1 

±  1 

±1 

p,q,r  rad. /sec. 

±1 

±  1 

±  1 

±1  j 

p,q,r  rad. /sec. ^ 

±4 

±  4 

±4 

±4 

/?  rad. 

±•5 

±•5 

— 

±•5 

X 

£  rad. 

±2 

±  2 

-- 

±2 

grad. 

±2 

±2 

— 

±2 

8  rad. 

±•5 

±*  5 

-- 

-- 

* 

8  rad. 

« 

±1 

±  1 

-- 

-- 

8  rad. 

±.125 

±.125 

•  « 

% 

±2 

±  2 

±  2 

±2 

±2 

±  2 

±  2 

±2 

0C  rad. 

±•5 

±5 

±  -5 

±  -5 

9C  rad. 

±•5 

±•5 

±•5 

±•5 

tic  rad. 

±•5 

±•5 

±.5 

±•5 

90  rad. 

±TT 

±  w 

±  t r 

±ir 

9  rad. 

0 

XX 

±4 

±  4 

±  4 

±4 

9  rad. 

0 

±4 

±  4 

±  4 

±4 

Q,~  rad. 

ks.T5 

±1 

±  1 

±  1 

±1 

9 S  7t-  rad. 

XX  *  ^ 

±4 

±  b 

±  4 

±4 

7r  rad. 

s  -75> _ 

±4 

±  4 

±  4 

±4 

TABLE  IX  (Continued) 


Rotor  Systems 

CH“3^ 

d 

f 

g 

(C ip/o-  )s 

±•5 

±•5 

±•5 

±•5 

(<V<r  ) s 

±  *5 

±•5 

±•5 

±•5 

(c y/a  )s 

±•5 

±•5 

±•5 

±•5 

(C1/°’)S 

±•0312 

±  .0312 

±•25 

±.0312 

)s 

±.0312 

±  .0312 

±  .25 

±.0312 

(<V  ^  >s 

±.0312 

± .0312 

±  .25 

±•25 

(Cn/<r  )f 

±.0312 

±.0312 

±.25 

±.0312 

(C j/er  )f 

±.0312 

±.0312 

±.25 

±.25 

)f 

±  .0312 

±.0312 

±.25 

±.25 

(Cg/o-  )e 

±.0312 

±.0312 

±.25 

±.25 

qW± 

±.025 

±  .025 

±•4 

±•4 

^w2 

±  .00625 

±  .00625 

±.l 

±.l 

± .003125 

±  .003125 

±•05 

±•05 

X 

,W1 

±.05 

±  .05 

±  .8 

±.8 

^2 

±.0125 

±  .0125 

±.2 

±.2 

X 

qw^ 

± .00625 

±  .00625 

±  .1 

±*  1 

\ 

±  .40 

±  .40 

±  6.40 

±6.40 

±  .10 

±  .10 

±  1.60 

±  1.60 

s 

±.05 

±  .05 

±  .80 

±.80 

qVl 

±.025 

±  .025 

±  .4oo 

±.400 

Qvg 

±  .003125 

± .003125 

±.050000 

±.050000 

X 

qvi 

±.050 

±.050 

±.800 

±.800  ! 

X 

qV2 

±  .00625 

±  .00625 

± .10000 

±.  10000 
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TABLE  IX  (Continued) 


Rotor  System 

CH-3C,a,b,c,e 

f 

g 

V1 

±.4o 

±  1.60 

±  6.40 

±6.40 

qv2 

±.050 

±.200 

±.800 

±  .800 

±.2 

±.2 

±.2 

±.2 

4, 

±8 

±8 

±8 

±8 

* 

±8 

±  8 

±8 

±8 
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Figure  79-  Hybrid  Computer  Flow  Diagram 
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Figure  80.  Continued 
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APPENDIX  III 


PROGRAM  EQUATIONS 

This  appendix  presents  the  equations  used  in  this  study.  The  fuselage 
equations  are  given  below  in  nondimens ional  form: 


X 

4* 

II 

- ^ 

<r  U 

_  1 

5*  s 

-  Kx 

K, 

On 

%> 

_  1 
M4. 

+  ip« 

P^f 

-  *>w,# 

V 

K* 

K, 

(£2 
v  0* 

+  £ 
k, 

& 

+ 

iW 

/«.*  = 

cjfUH 

-p/<* 

+ 

£a) 

<r 

+ 

K-% 

K, 

?*< 

Ki  L**  — 1 

k71v^ 

• 

p  = 

_ K59 

Ke 

tr  + 

K-,  (  Cy)  . 
Ke  vo-4 

k6 

V  <r 

!)f  + 

K*V>  <rk 

+X( 

K\x 

Ka 

jC*+ 

% 

=  —  r  p  —  4.  Kzo  (  ,  _L  (L.m) 

K,3l<r4  kTj  V  5=  /V  * 

^  K,‘  ^  K'»  fc  K.»  V  *»/&% 

-tKrV-^ 
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where  the  constants  (K)  are  defined  in  Table  X  and  the  axis  systems  are 
shown  in  Figure  8l. 


The  fuselage  angle  of  attack,  sideslip,  and  velocity  are  given  by: 


Fuselage  attitude  and  rates  are  defined  by  the  Euler  angles. 

fa  =  p  ♦  %  9< 

\y  s  r  cet>&4  ♦  ^ f 

^  cos, 
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o 
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The  fuselage  and  rotor  velocities  and  accelerations  are  related  by  means 
of  the  following  equations: 


■ 

\ 

=■ 

-  {yU*,  -  */a.(K,  ij 

f 

V*, 

= 

5 

A?y/ 

= 

-^YF  ^  -Ac  (  ^7  ~  K}0  l~) 

i 

a 

i 

&, 

- 

1+  + 

1 

Qy, 

— 

f/jlc 

i 

- 

jfcfr  T*  -  rJ 

] 

X 

r 

_  JJ.C  ^^.*F  S^iF^s  ~  K/ ~ 

1 

1 

; 

try, 

7b  C^i>"r  ■+  ib  ^  - 
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/irii  =  ~sr0  +  ^r/a0(Kie  +  k7Ij)J 

&*>  = 

Qy,^  =  */s£ 


The  rotor  flapping  and  lead  acceleration  are  given  by  the  following  equa¬ 
tions: 
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The  Z'  and  f'  terms  are  velocity  and  acceleration  components,  respectively, 
and  are  defined  by  the  following  equations: 

•rl  _  _  _  _ 

A  =  —  £di  +  a-s  ~^c„ai  +w,(a»  +ebjt^{ 

~fi  “  £^2  +  ^4^  +  dj  +  2.  §  d-2  i  |3[  b+  +  2  ft  hj 

"f«  =  b3  -  zfi  [at  +  ^3  (ft  4  “*,')]  -  30.4. 

a. ,  =  cuyicos  Yfi -Co,,  $IH 

0-t  =  COx,  cosYk  +  ujy,  SUJ  % 
a3  =  OJy,  COS  Yr  -  <£)„,  SIM  Tr 
0-4  -  COxi  COS  Y>  &yi  SIM  Yix 

Q-S  S  &oy,  COS  Yr  ~  /lTe>i,  SIM  Yr 

d6  =  faxt  cos  Y*  +  ‘trcy>  sim  Yu 

a7  =  frey,  COS  %  ~  trcx,  SIM  % 

0-&  =  iTox,  COS  Yk  +  iTey,  Siu  Yk 

b,  —  Si  y  co*, 

bx  *  TL  -v  Co*%  +  £ 

bs  —  Si  +■  CO*( 
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d,  —  ^mx 

da  *  £3x,  snY*)  +*  ^u,  c.°*  fn  +  V©^(SihYV 

-  £5*,[v:e^l  &«Yr  '  ^x.S'MYr  +e(25l+  co*,)] 

da  =  2  ?i[63Xt  co  *>  Y«  +co^(&'nYr]  —  Cu^coiYa  *»■  ajK|  s » w  Yr 

Blade  flapwise,  chordwis^  and  torsional  modal  amplitudes  accelerations 
are  given  in  the  expressions  below  where  the  C  coefficients  are  defined 
in  Table  XI. 

Flapwise: 


%Wi  ClL  ^((  ”  5AY5(Lnila]] 

io(^~^\c  +"f»>  (C4.-C7)tjli 

+  ^  ^  (C”  ~  C”'S*l3  +£  W^e  6!  -  t)  Cy„? 

+  4(U.)LjF  ,  siV*  c-kfJ 

+  t?,PMs fPJlcSl)T  -Ci5C!i;,f] 

p-l 

l 

-  Wft  51  'jvp  fV»3  C»oL  +  M4.C„l 

+  ?2  sin  ft]  CiliJ  +  mt 


■f2o  —  S  (  -  2  05^,31) 
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%'*l'  - 


'■  y  * 


Aiv 


4  =  /#'£*- 2S..JI  -2SJi -S*-/3l 

?s  —  b5  -  £0,,)]  -/3  a*  +  / 

—  — •  1  **  1 

'  z.  **  “f»  +  P  +fiSL 

mla  =  cos  A  -  e^sivj  ft 
mi*  =  Sin  ft.  4  ft,  cos  A 
MJ  =  s/yvft 

m*  =  X2  cos  ft  +  ^  SIN  ft 

m.  =  i/3#  (ii  +  «,,  +  l)  +  4  * 

w*  *  z©4  +  2(s(-a  +  arv*ih 

Chordwise: 


N 


C,4f  ^  |  /v-pa LS*YSa  m.a  4-  Sa^s*  mlctjj 

p  l  ^  fi'p1  t^o  (^n  ~  Cts)^?1  *'Cm  0-'b  “ 


C*°V 


**.  • 
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I 


m».Cirp  -m,Cl6p  -[^cosA  sim0.]cs3(, 

+  Or  *  £>  Mld  +  (£<>“  eSl)C2J? 

+  (42l  -Hl)Ci+p]  +  rn7  w  Wsrpj1 

-2(o.o*5)  (Zv?  j„r 

Torsion: 

**I  C  *vl  “  t&f"1  ^ 

—  + $  )  ^/^^».cMlo  +  C37 yi1 


1 

1 

1 


Terms  common  to  the  above  equations  are  given  by  the  following  expressions: 


?3  =  -2p>  (Ji  4  u3*,  +  I)  -(3  (fi  4-  Ati  4- 

-f*0  ~  [-A  +  60*,  4*  £  ] 

rtt*-  -  Z  sin  A  cos  00 


% 


1 


R 


it! 
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>t*v 


Ml  t  2c#,,  4  £ )  &  +*  ( SI  +  £3*,  4*  £  )  “I 

m9  «  l  -  2  sikj*^ 

?2?  *  f  il  +  4-  §  J  -  I 

0,4  +  (51  ♦  l)[-j M5(A+<f)  t/O-m*)] 

+  0*1  A  (  I  +  Wft) 


The  lag  damper  moment,  including  both  viscous  damper  and  a  spring  force, 
is  given  by  the  following  expression: 

Ml„  =  [c.3,  a  +  k.si*J  a]  $  ( /wi)f,0  aw. 

u  J  l«l  w 

-[k.co&A  -  C„4  sin  a)  t  (£*),.. 
-C0[<.os$,  X  (^VP)F.0 
-  SIN  6,^  (<Twt-)p=o  'fnj  ~  Co  £  -  Ko  & 


The  blade  root  end  shear  forces  include  aerodynamic  as  well  as  dynamic 
terms  as  follows: 


A*5a 


N 

-£1*7 

as  1 


Aw*.ai  i-  v6„  a,  -5b,1  *  ^  -  M8  (fc.)e  { 

**  *  2  •««  u  , 

-A®  1  -  -  St, 
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+  2.(  *a. 

-b,  fi7 


+  ^  sj)  ^  c 

'  4-  Tw/  '-lot 

+  s6)i§  c  4 

?*i  o^r  u  \sj> 


K*  >- 
^irp 


^  ^  ^7,  ^v? 


't  ^  ^  ^A**«-  “  Mb  (v.»,  -  e[a, 

~ <3-1  (SL  +  b,-)!  +  d,  +  5*j)  ~^B(^c«)ft[d*  +  co,(a, 

+/]  -  co4  4  t%.C„i  *  *,  t  4  .  C* 

C*l  t*|  0w4 

I  V>  c«r  -‘•fj.jC.p 

-^|^S,M.  -M. [a, 
-Vo,,!,  tw.,0,  +  eb,  ♦  §„j  -  Me(rt«)ft  p£ 
-l£ o.t  +  b,  -ib,(flj -t* S)  -  £  b,*} 

—  COi  6,  £  <JypC|Cf  +(S7  *  2**!S4) 

+  *b,  (sl+  ^  St)  Z  <*v? CiS?  +S.N  4  X  t«i  c ,oi 

T  1  l-l 
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-(2b, s,  -<>,) 

-  2b|(s>,  ^4s«)  £  a  w.  C.,oi. 

Ul 

S»  =  —  P  cos  A  +  J  St Kt  Qt 

x  v 

*  &  COS  0c  +  ft  SIN  Oo 

Si  =  — ^3cos 6c  4  5  SIN  & 

S4  *  £cos4  +/3S1vi4 
S«  =  (1  t  C1()  sm^# 

S4  •  (jl  +  C|.)  cos^ 

St  =  2  4  s«M<& 

«  2  4  cos  4 

The  rotor  blade  segment  forces  are  represented  as: 


^A*5a  Z  ^rnf  U(UTCJt  +  U?CD) 

^A'fSa,  ^  ^  ^m7  ^  ^  *  UT  C  6^ 

M«aAa.  =  CmeM)Ua  +  ^*Sa 
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K 


-  "iv Ut**° 

=  >J/4C  C *»*■  %  *  &)  Wt»-  <  0 


Blade  segment  velocities  are  given  by  the  following  equations: 


U, 


=•  "Ao  ( I  —  \  )  4-.  /3  (a76  -v-  a.p)  +  a,(e  +rn) 

-r«.  (ax  £  +  p)  - 

4  |r/m  V7*  4  Vv-cc)  “  ;lfi  ^  ^wc<t  U] 


+  SiN  A 


?=/ 


ae(^i^T,a  Vp  +  >»Mcjj  '')v K%<-C%e° 

uT  =  <*-i  (»-£■)  -  as£  +  (e  +?*.)[£ -|3 {a,* 51 4)] 

2  r"  2L 

.+  (5  *  rOk,  +  CCS  6.  ( try,  L-r  -  a  J  JJ  x 

?=i  SJI  ; 

i 

_«  1  3  ^  , 

tA3%u.c.])  ^wt  ~  a5 


U  —  H'p  S*N  4*  \X ^  COS 
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where  the  blade  segment  angle  of  attack  is  given  by: 

o<r  =  &.  +• 


*  -  ™'(&) 

and  $0  is  the  actual  geometric  pitch  of  the  blade  at  the  three  quarter 
radial  station,  r/R  =  0.75* 
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9, 


(9i\s  -tm  S,j>  +Z(M?.0  ^-t] 

+  J:  ^  *J  v  ^  m  +  ^  (^p*o  i  «t  ] 


4 


TAW 


T«i 


The  inelastic  blade  pitch  is  given  by  the  following  equation: 


"[A.sCOSYr  t  B,4SIN  Y^]  +.15  ft 

C^O.is  “  [. Ats  +  B»s5l]  coi yR  t^Ais 5i 

-C>»b]  s>w  Y* 

^).75  ~  4[(A|S&  -  -  A,*  - 

*[(2  A»b  +  B»sJl)Jl-Bis  +Ai*jl]s»N  Ya 


The  pitch  of  the  blade  along  the  blade  radius  is  given  by: 


56  =  (fL  +  €  -  .75)  e, 


The  various  gravity  vector  terms  are  defined  as  follows: 
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Rotor  shaft  forces  and  moments  are  given  by  the  following  equations: 


—  y  c 
|°  bcR  4"  b*s 


SB 

4 


These  forces  are  then  transformed  into  fuselage  forces  and  moments. 


(H  =  (H 

(H= 


r  ^ 


Vt 


*■'  1  ■■■— I1>WW  J 
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(H  -  (H  -  ( 

(ty  -  (r), 

(H  *  (?').  *  (Hi. 


The  nondimens ional  rotor  speed  is  calculated  using  the  following  equation: 


il(nT) 


.TUnT-T)  +Sl„Ta(nT) 


where  the  acceleration  of  the  rotor  is  calculated  using  the  following 
equation: 

9. 


LIhu**‘"-R  be"  Mb] 


mc 


<3. 


— 


Prior  to  the  simulated  loss  of  the  rotor  driving  torque  the  following 
equation  is  used: 


At  the  time  of  the  simulated  power  failure,  one  of  the  following  two 
equations  is  used: 


220 


Full  power  loss: 


(Cq/<r)e ( nT)  -  (TisJ) (c*}/<r)e ( nT -t) 


Part  power  loss: 


Cc,/<r)e(nT)  =  ( '^=J)  (c9/(r)e  (nT-T)  +  (fjGC,/r). 


The  rotor  blade  azimuth  station  part  power  loss  is  calculated  using  the 
following  equation: 


VkCnT)  =  +SI.T1 t  (nT) 


The  local  blade  bending  moments  were  calculated  using  the  following  equa¬ 
tions: 
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M.  = 


Other  quantities  calculated  were  the  chordwise  bending  moment  in  the  vicin¬ 
ity  of  the  lag  damper. 


and  the  deflection  of  the  tip  of  any  blade  within  5  degrees  azimuth  posi¬ 
tion  of  the  tail: 


+  wC), 


Tail  Rotor  Equations 

The  tail  rotor  thrust  is  calculated  using  the  following  Bailey  theory  equa¬ 
tions,  which  include  the  effect  of  delta  three: 


~  Y'  “ 


~  ^57 


Two  additional  equations  were  used  to  calculate  the  change  in  altitude  fol¬ 
lowing  a  loss  in  power  and  normal  load  factor,  and  are  given  by  the  follow¬ 
ing  equations,  respectively: 
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A  2  =  &  J  '&***  $4.  '***'  4+  "^^94  Ot*>  ${]co*6^(jjt 

and 


^  -^2.# 

7 


NONARTICULATED  ROTOR  CONVERSION 


The  above  equations  were  modified  for  the  hingeless  rotor  and  the  floating 
hub  portion  of  this  study. 

Hinge less  Rotor  Modification 

The  lead-lag  and  flapping  equations  were  eliminated,  and  flapping  angle  was 
replaced  by  a  fixed  coning  angle. 


-/■O 


X  XX 

/3-  P  *0 
P'  Pb 


The  rotor  acceleration  was  calculated  using  the  following  equation, 


A 

Jh 


[l»  F'b  (I*+ 
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^X2)J 


and  the  blade  segment  velocities  were  modified  to  result  in  the  following 
equations: 


+(*+**  )[*>,  7s  fa  +  T^)J 

+  l*)  5a~-J>64  6, 

Blade  flapwise,  chordwise,  and  torsional  modal  amplitude  accelerations  were 
modified  to  include  the  following  effect  of  Coriolis  acceleration  and  struc¬ 
tural  damping: 

To  the  right  hand  side  of  the  flapwise  equation  the  following  term 
was  added: 


cu  C~  cn-i(ix  Ja  #.)] 

The  quantities 

C/4*  L  ^ 2 -f'1  A  ^  <£.„■'  )~  ^47 f(46  Cc~<- 

+  £  9'  )]-  0.65 

were  added  to  the  chordwise  equation,  and  to  the  Sy^  equation,  the 
following  term  was  added: 

C  t0  1 1  £>t  (z. .  SCp  ^  Bo 

The  moment  at  the  root  of  the  blade  was  calculated  using  the  following 
internal  moment  equations: 
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•M 

« 

—  r  —? 

M  ^  3  =  O-  75  S,  &o]  ^  C$i  ^  ^  -^b 

~  E***  @0  ~ons&t  Cr*+  &JJ  ^  Csisi. 

-  -  C^~*  &6  +  O.ISG,  6[J  2  C-sisC  ft! sotrx 

-  ^  £5/40 

so  that  the  rotor  moments  can  now  be  written  as: 

fel- 

(CA' 

(A)’  - 

Floating  Hub  Modifications 

The  floating  hub  rotor  system  was  represented  by  two  orthogonal  teetering 
rotors. 

The  flapping  equation  was  changed  from  the  previous  equation  to  the  follow¬ 
ing: 
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Before  integration,  the  flapping  accelerations  for  opposing  blades  were 
made  to  conform  to  the  teetering  condition: 

**  .  v  X*  At  v 

/Sj  £y)  &  (2l±J3°1L 

=-is\  V 

J&X  +  96')  -  A  fct±  (iL ±2J0‘) 

2.  " 

At  ('*+Z  10*)  -  -f3t  y  cfo  •) 

The  flapping  velocity  and  angle  for  the  floating  hub  then  becomes, 


j|  ^ 

/St  ('f)  •  M  (’v)cU 

A  ('Y*  iso*)  --  -  fA  e+)dt-A  (v) 

A  (V+  90 •)  --  j  A  (V'  +  9o‘)cbt 

A.  (Y +216*)  --  ~f/3t  (V  +  90’) - -/SjV' +9*") 

fiCP)  -|/3 '(*)*&+& 

(i  (y  +/SO*)  =  -P(*)  +Pa 

(i(y+  90‘)  *  JA  (V  +  90‘)ctf 

fi  /V  +110°)  =  ’ft  (Y  +  90*)  +  (2t 
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where  !•  a  built-in  coning  angle. 

Blade  segment  velocities  were  expressed  as  follows: 


K  ( >-  T)‘hfi(<*ts)  +  (i-h'Y^)Ectl  -/sj 

Sa  ^  (■&  7a  Cj q  ~] 

+  L  ^a  -^4-a  'hJ>i.4  J 1  <£t 

-^r  *  )  C 4  -p(at  +TL  %  )7 

V-  Z*-^3q  ~]  Cc-a.  Q.~  [j>Sa  <9e 

For  the  floating  hub  configuration  Vx oh out  the  elastic  degrees  of  freedom 
the  rotor  moments  become: 

-  -jftx 

(c>/r\  -  ° 
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TABLE  XI 


C  CONSTANTS 


The  blade  constants  appearing  in  the  various  equations  of  motion  are  de 
fined  below: 
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